A new approach for doing theoretical and numeric work with 
Lemaitre— Tolman— Bondi dust models. 
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We introduce quasi-local integral scalar variables for the study of spherically symmetric Lemaitre- 
Tolman-Bondi (LTB) dust models. Besides providing a covariant, and theoretically appealing, in- 
terpretation for the parameters of these models, these variables allow us to study their dynamics 
(in their full generality) by means of fluid flow evolution equations that can be handled with simple 
numerical techniques and has a significant potential for astrophysical and cosmological applications. 
These evolution equations can also be understood in the framework of a gauge invariant and covari- 
ant formalism of spherical non-linear perturbations on a FLRW background. The covariant time 
splitting associated the new variables leads, in a natural way, to rephrase the known analytic solu- 
tions within an initial value framework in which covariant scalars are given by simple scaling laws. 
By using this re-parametrization of the analytic solutions, we re-examine and provide an alterna- 
tive outlook to various theoretical issues already treated in the literature: regularity conditions, an 
Omega parameter, as well as the fitting of a given LTB model to radial profiles of density or velocity 
at different cosmic times. Other theoretical issues and numeric applications will be examined in 
separate articles. 

PACS numbers: 98.80.-k, 04.20.-q, 95.36. +x, 95.35. +d 



I. INTRODUCTION. 

Spherically symmetric LTB dust models are among 
the best known exact solutions of Einstein's equations 
Since they provide access to non-linear effects as- 
sociated with inhomogeneous sources by means of ana- 
lytic and/or tractable numeric solutions, they have re- 
ceived widespread attention in the literature (see 0, 01 
for comprehensive reviews). These exact solutions are 
used mostly to study cosmological inhomogeneities, as 
for example the series of articles in the references 0- 
[loj , but have also been employed to examine a variety of 
issues, for example as useful toy models to describe grav- 
i_t_ational collapse under a classical |3, IB| and quantum 
aroach. LTB models have also been widely used 




llHlSj in the context of the ongoing widespread theo- 
retical and empiric effort to explore the possibility that 
cosmic observations could be explained without resorting 
to dark energy, but by taking into consideration the fully 
non-linear effects of inhomogeneities in the context of 
General Relativity (see [l^ for a review). These models 
are also helpful as theoretical tools to probe various aver- 
agirig techniques applied to cosmological inhomogeneities 
|14l [la |2(]| - |22| ■ in particular, the scalar av erag ing formal- 
ism developed by Buchert and coworkers 23|. One finds 
in the extensive literature on these models a preferred 
set of free functions and analytic solutions (implicit and 
parametric) that has become a sort of standard 0, H, [13] 
(even in numeric work on these models [T^ - [l8| b 

Given the success of these models and the fact that 
their standard parametrization is adequate and works in 
practice, then it is legitimate to ask for the justification 
of an article proposing and discussing alternative vari- 
ables. The answer to this question is simple: alternative 



variables may provide or motivate either new theoretical 
and empiric developments or allow us to grasp known 
results under a different perspective, thus illuminating 
possible unsuspected connections with other models or 
theoretical issues. In this context, we propose an alter- 
native description of LTB dust models that is based on 
a set of covariant scalar variables defined by applying 
to arbitrary scalars the same type of integral construct 
that yields the Misner-Sharp quasi-local mass function 
|26l - l28l |. Hence, we find it natural to call these variables 
quasi-local scalars. 

The quasi-local scalars mentioned above have already 
been used for a dynamical system study of dust LTB 
models [1^ and in the application of Buchert's averaging 
formalism to LTB models [21, 22], as well as in a numeric 
study of inhomogeneous dark energy sources associated 
with an LTB metric, but with an ener gy- momentum ten- 
sor that is more general than dust [30l l3l| . The pure dust 
case merits the separate study that we provide here, not 
only because dust is an adequate source to model cold 
dark matter in cosmological scales (and a much less ar- 
tificial source for an LTB metric than the fluid tensor in 
;')(). _31]), but because dust sources allow for a qualitative 
analytic framework that complements numeric work and 
yields very useful theoretical and practical information. 

The potential of the quasi-local variables for ana- 
lytic and qualitative work on LTB models follows read- 
ily from parametrizing their known analytic expressions 
in terms of these variables and their fluctuations. The 
new variables lead in a natural way to an initial value 
parametrization of LTB models, which in turn leads to 
the introduction of a FLRW-like metric and scaling laws 
for density, Hubble expansion, spatial curvature, shear 
and all other covariant quantities of physical or geomet- 
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ric relevance. Various elements of this parametrization 
have been already considered in the literature [l3l - [T7| . 
but only as a useful coordinate ansatzes that are jus- 
tified because the involved expressions resemble FLRW 
parameters. In this article we provide an appealing the- 
oretical and covariant context for these expressions that 
could allow for new empiric or numeric results and/or a 
better understanding of existing work. 

The quasi-local scalar representation applied the LTB 
dust case leads to a complete description of the dynamics 
in the context of "fluid flow" (or "l-(-3") evolution equa- 
tions that are fully general, but can be handled by simple 
numerical methods, and thus have a significant potential 
for empiric modeling of cosmological dust sources. This 
scalar representation also leads to a description of LTB 
dust models in the framework of spherical, non-linear, 
gauge invariant and covariant (GIC) perturbations over 
a FLRW background (in the context of the fluid flow 
perturbation formalism derived by Ellis and Bruni ^32| - 
|34|). The resulting perturbation framework is also con- 
sistent with the traditional perturbation approach "s?], 
as the evolution equations of quasi-local scalars and their 
fluctuations generalize the evolution equations of linear 
perturbations of spherical dust sources in the comoving 
gauge [38]. 

The plan of the article is described below. We sum- 
marize in section II the metric, field equations, analytic 
solutions and basic regularity of LTB models in the con- 
ventional variables. In section III we introduce the rep- 
resentation of local covariant scalars and their evolution 
equations in the context of the covariant time slicing asso- 
ciated with a fluid flow approach to the dynamics of the 
models. The quasi-local scalars and their fluctuations 
are introduced and discussed in section IV, while section 
V deals with some theoretical and practical implications 
of these scalars. In particular, the quasi-local scalars 
lead naturally to an initial value approach in which ini- 
tial value functions are given by the new variables eval- 
uated at a fiducial (initial) hypersurface defined by the 
I-I-3 time slicing. This approach leads to a FLRW-like 
form for the LTB metric, which has been suggested previ- 
ously (see but as an ansatz without a theoretical 
context. The initial value approach also leads to simple 
scaling laws for all scalars and an intuitive gauge for the 
radial coordinate. Hence, in section VI we rewrite the 
analytic solutions of the models in terms of the above 
mentioned initial value parametrization, which leads in 
section VII to an elegant characterization of the coor- 
dinate locus of the curvature singularities that is sim- 
pler and more intuitive than that using conventional vari- 
ables. In section VIII we rewrite the Hellaby-Lake con- 
ditions [U, [3^ to avoid shell-crossings in terms of the 
initial value functions, while in section IX we rewrite the 
analytic solutions in the same initial value approach, but 
using a suitably defined Omega parameter and the quasi- 
local expansion factor as initial value functions (which 
can be used also to parametrize the Hellaby-Lake condi- 
tions). 



In section X we address (as an application) an inter- 
esting theoretical issue raised by Hellaby and Krasinski 
(see 0, [13, Ell), namely: the possibility of mapping by 
an LTB model of a given density or velocity radial pro- 
file at one cosmic time to another density or velocity 
profile at a latter time. We show how the initial value 
parametrization of the analytic solutions greatly simpli- 
fies the understanding and handling of this problem. In 
section XI we derive three equivalent systems of fluid 
fiow evolution equations for the quasi-local scalars and 
their fluctuations. Each of these systems provides a full 
description of the dynamics of the models, and all are 
effectively equivalent to systems of ordinary differential 
equations subjected to constraints, hence they are tech- 
nically more accessible than (and as fully general as) the 
system that arises for the local scalars in the H-3 fluid 
fiow framework. We explain how initial conditions can be 
set up and argue that these equations have a significant 
potential for numeric computation of quantities of ob- 
servational interest (null geodesies, red shift factor, area 
and luminosity distances, etc). In section XII we show 
how the quasi-local scalars, their fluctuations and their 
evolution equations can be understood in the framework 
of a gauge invariant and covariant formalism of spherical 
non-linear perturbations. 

In section XIII we derive a fourth system of fluid flow 
evolution equations, but now for the Omega and Hub- 
ble parameters introduced in section IX. These evolution 
equations lead to a dynamical system (like the one in 
[29|). We comment on the relation between the perturba- 
tion formalism of section XII and theseparameters (and 
others proposed in the literature [H, [l3l), thus provid- 
ing a theoretical context for quantities given as coordi- 
nate ansatzes on the basis of their resemblance to FLRW 
Omega and Hubble parameters. We discuss in section 
XIV various aspects of special LTB conflgurations, such 
as: closed models, models with simultaneous big-bang or 
collapsing singularities, simultaneous maximal expansion 
and a mixed elliptic/hyperbolic configuration. Our con- 
clusions and suggestions for future work are summarized 
in section XV. In Appendix A we discuss various regular- 
ity issues of the quasi-local scalars and their fiuctuations, 
while Appendix B discusses the relation between quasi- 
local spatial curvature (which determines the kinematic 
class: parabolic, hyperbolic, elliptic) and local spatial 
curvature. 



II. LTB DUST MODELS IN THE 
CONVENTIONAL VARIABLES. 

Lemaitre-Tolman-Bondi (LTB) dust models [ll, are 
the spherically symmetric solutions of Einstein's equa- 
tions characterized by the LTB line element and the 
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energy-momentum tensor 



= -c^df + dr^ + i?2(d6'2 + sin^ Odip^), (1) 

1 + 

T"^ = pc^u''u\ W^S^, (2) 

where R = R{t,r), R' = dR/dr, E = E{r) and 
p = p{t, r) is the rest-mass density. The field equations 
reduce to 



• , 2M 
2M' = npR^R', 



(3) 
(4) 



where k = SttG/c^, M = M(r) and 7? u'^Vai? 
dR/d{ct). 



A. Kinematic equivalence class and analytic 
solutions. 

It is common usage in the literature (see 0, S [13]) 
to classify the solutions of ([3]) in "kinematic equivalence 
classes" given by the sign of which determines the 
existence of a zero of R? . Since E = E{r), the sign of 
this function can be, either the same in the full range of 
r, in which case we have LTB models of a given kine- 
matic class, or it can change sign in specific ranges of r, 
defining LTB models with regions of different kinematic 
class (see [24| ) . The solutions of the Friedman-like field 
equation ^ for each kinematic class take the following 
well known parametric form: 

Parabolic models or regions: E = 0. 

c(i-i,,) = ^?7', R^{2My/^rj-\ (5) 
Hyperbolic models or regions: E > 0. 



E = — (cosh ?7 - 1) , c{t-t^^) ^ (smh ?? - ry) 



Elliptic models or regions: E <{). 



M , 

E = — (1 ~ cos ?7) , 



M 



|£:p/2 



(?/ - sin 77) , 
. (7) 

where i^b = ^bbir) is called "big bang time", as it marks 
the coordinate locus of the central expanding curvature 
singularity: R{t,r) = for r > 0. It emerges as an "in- 
tegration constant" in the integration of ([3]). Notice that 
the locus of central curvature singularity is distinct from 
that of the center of symmetry R{t, 0) = (see Appendix 
lA Besides the kinematic class, LTB models that ad- 
mit (at least) one symmetry center can be classified as 
"open" or "closed" , respectively corresponding to the hy- 
persurfaces of constant t being topologically equivalent to 
or §3 (see Appendix [X3)) . 



In order to deal with a given specific model by means 
of ©-([T]), we need to prescribe the three "conventional" 
free functions 



M(r), E{r), t,,{r). 



(8) 



However, the metric ([1} is invariant under rescalings of 
the radial coordinate r — r{f), hence it is always possible 
to eliminate one of the free functions in by a suitable 
choice of the radial coordinate, thus leaving only two ba- 
sic irreducible free functions. 

Given a choice of free functions ([H]), the solutions ([5])- 
([7]) can be used to find the remaining relevant quantities 
that may be required for a specific problem, a procedure 
that has been used abundantly in the literature [3, Q . 



B. Regular LTB models 

We will assume henceforth the existence of (at least) 
one symmetry center marked by r = (see Appendix 
lA and will denote by "regular LTB model" any LTB 
configuration for which the condition 



sign(i?') = sign(il/') = sign(Vl + £^), 



(9) 



holds for all t and all r not marking a symmetry center 
(since M'(0) 0, E{0) = and R'jt 0) = 1 [mHI). 
Conditions © imply that if R', M', \/l + E have a zero 
under regular conditions, it must be a common zero of 
the same order at some worldline r = rtv > 0. Violation 
of jSl), I.e. R' = for M' ^ 0, ^\ + E ^ 0, resuhs in a 
shell crossing singularity for which p — > cx) occurs with 
i? > 0. Condition ([9]) together with (|4]) also implies 
that p > holds everywhere and that it is bounded 
everywhere except at the coordinate locus of a central 
singularity (i.e. points for which R(t,r) — holds 
that are not the worldline of a symmetry center). The 
restrictions on the free parameters ([8]) that guarantee 
the fulfillment of ^ are the well known Hellaby-Lake 
(necessary and sufficient) conditions [H,!!!,!!!] given by 



(6) Parabolic and hyperbolic models or regions: 



i?' > ^ {Af ' > 0, E' > 0, C < 0} , (10) 
Elliptic models or regions: 



± i?' > ^ 



±M' > 0, < 0, 



± 



2 E 



27rM 



>0, 
(11) 



where only expanding configurations are considered in 
([TU| and the ± sign in PTI) accounts for the fact that R' < 
occurs in elliptic models whose that admit a second 
symmetry center (see Appendix IA3|) . The equal sign 
holds only at symmetry centers and at values of r where 
R' = 0. 



4 



III. COVARIANT OBJECTS AND THEIR "1+3" 
EVOLUTION EQUATIONS. 

The normal geodesic 4-velocity in (U)-© defines a 
natural time slicing in which the space slices are the 
hypersurfaces ^T[i], orthogonal to and marked by 
arbitrary constant values of t[4^. The metric of the 
'^T[t\ is simply hab = 9ab + UaUh = gtjSldl (or ^ with 
dt = 0). The proper volume element and 3-dimensional 
Ricci scalar associated with these hypersurfaces are 

dVp = ^det{hab)d^x = == ^, (12) 

V 1 + ^ 



2{ER)' 
R^R' 



2E' 2E 



R'R i?2 



(13) 



Given this time splitting, the covariant objects in dust 
LTM models are p and ^TZ in ^ and (HS]), together with 
the shear and electric Weyl tensors, (Tab — ^ (aUb) — 
(e/3)/ia& and EabUcUdC^"'^ , where Va = h^Vb, and 
(jabcd ^j^g Weyl tensor, and the expansion scalar, 6, 
is given by 



2R 
~R 



R'' 



(14) 



Since LTB models (as all spherically symmetric space- 
times) are locally rotationally symmetric (LRS), all their 
covariant objects are either scalars or entirely charac- 
terized by covariant scalars [s^. Therefore, the spatial 
trace-less shear and electric Weyl tensors can be given 
in terms of single scalar functions, S, f , in a covariant 
manner: 



ab 



■~ab 



1 




ab 



6 



(15) 



where E"'' = ft,"^ - Sx^'x'' and x° = Vh^ is the unit 
vector orthogonal to m° and to the 2-spheres orbits of 
S0(3). 

A description of the dynamics of LTB models, which 
is an equivalent alternative to the field equations and 
their analytic solutions ©-([7]), follows by considering 
the following set of covariant scalars 



A = {p, e, 37^}, 



(17) 



which completely characterize these models. Hence, their 
dynamics is completely determined by suitable evolution 
equations and constraints for these scalars. Following 
the "fluid flow" or covariant "1+3" framework of Ehlers, 
EUis, Bruni, Dunsbury and van Ellst [US HI], and 
considering (|15p and (IT51) . the 1+3 evolution equations 



for the scalars p7)) are 



t = -— S + S2_£, 

3 

5.-|,E-3£(| + E 



(18a) 
(18b) 
(18c) 



(18d) 



while the spacelike and Hamiltonian constraints are 



e\' R' K , , R' 

S+- +3S-^0, -p'+£' + 3£-^0. 



3 =3^-^ + ^' 



(19) 
(20) 



Solving this system of partial differential equations de- 
termines the metric functions and remaining quantities 
though the definitions ©, (g]) and ([12])- (HH) that relate 
these functions to the scalars in ([TT]). However, it can 
be a difficult system to use in numerical work because 
time and radial derivatives cannot be (in general) de- 
coupled. In the following section we will introduce an 
alternative set of covariant scalars whose corresponding 
evolution equations (see section XI) are wholly equivalent 
to (fT8t-(l20l) but are easier to handle. 



IV. 



QUASI-LOCAL SCALAR FUNCTIONS AND 
FLUCTUATIONS. 



The function M that appears in ©-(g]) is for LTB 
models the Misner-Sharp quasi-local mass-energy func- 
tion, which is an important invariant in spherically sym- 
metric spacetimes [26l - t28| . It is basically the integral 
along a spherical comoving domain of the field equation 
(|4]), which will be well defined if we assume the existence 
of a symmetry center (at r — 0) and can be given as 



2M 



2G 



[ pTdVp^K [ pR^R'dx, (21) 

Jv(r) Jo 

(22) 



j' = VTTe 



i?2 + 1 



2M 



R 



1/2 



where the spherical comoving domain Vir) is defined fur- 
ther below and we have used the notation ...dx — 
Ji^~^ ...dx. It is evident from ([21]) and ([22]) that M is a 
proper volume integral "weighed" by the invariant scalar 
which generalizes the "7" factor of Special Relativity 
0,[2i|. The integral definition motivates the in- 
troduction of similar "weighed" proper volume integral 
functions which turn out to be very useful in the study 
of LTB models. In order to define such integrals, we re- 
mark that every regular '^T[i] contains comoving regions 
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(containing a symmetry center) that are diffeomorphic to 
the product manifold 



so that 



d{r) = {a;|0 < a; < r}, (23) 



where x = marks a symmetry center. Because of spher- 
ical symmetry every scalar function on a domain I?(r) is 
equivalent to a real valued function on the real interval 
z?(r). 



A. Quasi— local scalar variables. 



Consider now the following 

Definition 1: Quasi-local scalar functions Aq . Let t?(r) 
be a radial domain in a regular '^T[t]. For every scalar 
function A[to] : i?(r) R, there is a "quasi-local dual 
fimction" Aq[to] : '!9(r) M such that 



J^AJ'dVp _ J^'^^ A(to,x)R'^{to,x)R'{tQ,x)dx 



/^''J-dVp R^ito.x)R'{to,x)dx 

(24) 

where J- is given by p2|) . Since it is clear that i is an 
arbitrary but fixed parameter in these integrals, we will 
omit (unless it is necessary) to express explicitly the 
time dependence of scalars A and Aq . 

Comment 1 : Definition (fM)) is invariant under arbitrary 
rescalings of the radial coordinate r — r(f) that do not 
violate ([9]) and ()A7|) . In particular, the proper length 
^ = / y^grrdr Can be used as the integration variable at 
each '^T[t] (though it is not useful as a global integral 
variable because £ = £(t, r)). 

Comment 2 : In general, not all hypersurfaces '^T[t] in 
LTB models are fully regular, since the central singular- 
ities (expanding or collapsing) are not simultaneous and 
will intersect some of the ■^T[t]. For these '^T[t] the do- 
main (|23p must be suitably restricted and the integrals in 
([M)) must be treated as improper integrals. We discuss 
this issue in IA51 



Comment 3. The quasi-local functions Aq in (|24|) are 
not "averages" , as they do not comply with the prop- 
erties of average distributions of a continuous random 
variable. They can be recast as proper volume averages 
with "weight factor" J- if we define them as functionals. 
See [231] for discussion and clarification of this issue. 



B. Properties. 



The definition (|24|) leads to the definition of a "quasi- 
local volume" given by 

V,(r) = [ TdVp = ATT I R^R'dx ^^R^{r), (25) 

Jv(r) Jo 3 



Va' R~ 



v 



3i?' 
~R 



where we expressed ([T^ as = [ln(i?^i?')] ' and the com- 
mutation oi d/dt with the integrals in ([M)) . The quasi- 
local functions comply with the following properties: 



A'=iAqy^^[A-Aq] 



1 



A'Vqdx, 



Aq^iAqy^{A)q + {QA)q-eqAq, 



(27a) 

(27b) 
(27c) 



C. A representation of covariant quasi— local 
scalars. 

The introduced quasi-local variables provide a com- 
plete and very useful representation of covariant scalars 
that is alternative to ([T7)) . Applying the definition ([M)) 
to ([3]), (fT3)) and (|21l) we obtain the quasi-local duals of 
p, ^TZ and 6 



2M 



k - -— V"^ - — - Om 



where, to simplify notation, we have defined 



2m= |p, 



k EE 



6 ' 



e 



-kq. (28) 



(29) 



Notice that M, R, R = u°Vai? and J" = VTTE are 
invariants in spherical symmetry (26l - [28j . hence niq, kq 
and Hq are covariant quantities. As a consequence of 
([28l) , the scalars E and £ in (fTS]) and (fT6|) , associated with 
the shear and electric Weyl tensors, become expressible 
as deviations or fluctuations of the local scalars H, m (or 
QtP) with respect to their quasi-local duals: 



iU-Hq), £ 



(to 



(30) 



Since the shear and electric Weyl tensors are covariant 
objects, then these fluctuations are also covariant. 



D. Fluctuations of quasi— local duals. 

A convenient way to relate the local covariant scalars 
A in pT)) and their quasi-local duals Aq follows by intro- 
ducing the relative deviations (fluctuations) 



A-Aq 

A„ 



Aq/Aq 
3R'/R ~ A, 



\- [A'Vqdx, (31) 

qy) JO 



where we used (j27a|) and (j27b|) . Considering (pO| and 
(pri) . each scalar A in ^T7\ can be uniquely expressed in 
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terms of its associated Aq and S^"^^: 



m 



1 + ,5(") 



1~L — 1~Lq 



(32) 



A. A FLRW— like metric and scaling laws for 
covariant scalars. 



Consider the dimensionless scale factor 



(36) 



By applying the definition (j3ip to the Friedman-like 
equation for in we obtain the following important 
constraint 

2?^2j(«) = 2m, - A:, (5('=\ (33) 



while the relation between (5'™^ and (5'-'°-' and radial gra- 
dients of the conventional variables M and E 



M' 



3R' 



1 + (5('") 



'e 



m' 



5(k 



(34) 



follows readily from (gll) and ([5T|) . 

Since to, fc, "H and m,, kq, Hq are all covariant scalars, 
their j(^) are also covariant. In fact, these fluctuations 
can be understood in the context of a covariant gauge in- 
variant and non-linear per turbation formalism (see sec- 
tion XII and references [S^, HH ) . Hence, the set 

{Aq,6^^^} = {mq, kq, Tiq, 6^^\ S^''\ S^^^}, (35) 

constitutes a complete covariant scalar representation 
that is equivalent and alternative to (IT71) . 

There are important regularity issues associated with 
variables psp . which we discuss in detail in Appendix A. 
For example, any one of the (5'"^^ might diverge if 
under regular conditions (no violation of ^ and finite 
curvature scalars). However, as we show in Appendix 
lA 41 this does not lead to a curvature singularity, as the 
Riemann tensor frame components do not diverge (see 
Appendix lA 2[) . Also, the regular radial range of some 
of the hypersurfaces ^T[t] is necessarily restricted by a 
central singularity, which does not invalidate the regu- 
larity of these scalars away from the singular coordinate 
surface (see Appendix I A 5|) . While these issues must be 
taken into consideration, they do not affect the usage of 
the variables ([55]) . 



V. SOME THEORETICAL AND PRACTICAL 
IMPLICATIONS OF QUASI-LOCAL SCALARS. 

The quasi-local variables and their relative fiuctua- 
tions introduced in the previous section lead in a natural 
way to an initial value parametrization of LTB models, 
which can be very helpful for numeric and qualitative 
work. Various elements of this parametrization has been 
already considered in the literature [T^-[T7|. but only as 
a useful set of coordinate ansatzes whose justification is 
merely their resemblance with FLRW parameters. How- 
ever, these ansatzes can acquire a clear covariant meaning 
by its relation with quasi-local scalars. 



where Ri = R{ti,r) and t = ti marks a fiducial initial 
hypersurface ^% = '^T[ti\. We shall denote henceforth all 
quantities evaluated t = U by the subindex i . Bearing 
in mind ((36l) . we can re-write ((28l) as 



^^qi J kql 



Hi 



L2 

L2 



^TJTjq kq 



2to„ 



kqiL 



L3 



(37) 
(38) 



which, if we identify L with a position dependent FLRW 
scale factor, are formally identical to the scaling laws for 
density, spatial curvature and Hubble factor in FRLR 
dust universes. Applying this parametrization to the 
LTB metric ^ yields the FLRW-like line element 



flf dr^ 

1 kqi 



Rj (d02 ^ gjj^2 ^^^2^ 

(39) 



where the new dimensionless metric function F is 
^_ R'/R ^ , L'/L 



R'^|R^ 



1 



R'JR^' 



r. = 1. 



(40) 



The local density (to = np/S) and spatial curvature (fc — 
^TZ/Q) satisfy the following scaling laws in terms of L, F 
and initial value functions 



+ ^(™)l ^ ^ 



k -t- ^[I + ^W] 



3(l + <5f)) 



(41a) 



(41b) 



where we used (g]), and Comparing (gT]) with 

([57|. and bearing in mind (^5)) and ([5Tt . we obtain the 
following scaling laws for i^^™^ and 6'^'^'^ 



. ^(™) 



5 + *'" - 

2 



1 + 5, 



(m) 



F ' 
2/3 + 



{i/2)5. 



[k] 



Applying (I3ip to the Hamiltonian constraint 
the scaling law for (5'^^ 



(42a) 

(42b) 

(42c) 
yields 



2j(w) = 2to, - kq 5^=) _ 2to,, - kq, L 5(>''> 



'^f^qi kqi L 



(43) 
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Inserting (|43l) into (l30l) allows us to find scaling laws sim- 
ilar to the ones above for H — 30 and the scalars S, f , 
respectively associated by (fT5|) and ([T6| to the shear and 
electric Weyl tensor. Using (PT|) . we can express the scal- 
ing laws (l42t as relations between the gradients H'^, m'^ 
and k'^ and L'. Notice that it is always possible to use 
([55)1 and to eliminate any one of the three pairs 
Hq, (5'-^' or rriq, (j'™-* or kq and S^''^ in terms of the other 
two. 

The scahng laws (gH), (gH), and (gS]) depend on both 
L and F. While the effects of their dependence on L is 
easy to grasp because the qualitative time dependence 
of L (for a fixed r) can be appreciated directly from the 
Friedman-like equation ([38]), there is no simple way to 
guess the qualitative behavior of T (cither in the t or r 
directions). As a consequence, there is not much we can 
do with these scaling laws as long as we lack an expression 
for r in terms of L and the initial value functions (we 
obtain this expression in section VII). 



B. A radial coordinate gauge. 



The metric of LTB models in the form ([39)) is also in- 
variant under an arbitrary rescaling r = r(f). Since L 
and r are time dependent, it is natural when using ([39]) 
to identify this radial coordinate gauge freedom with the 
freedom to choose the initial value function Ri{r). How- 
ever, because of \AQ\ and ([A7p . the choice for a function 
Ri is not completely arbitrary: it depends on the topol- 
ogy of the space shces '^T[t] (see [2^ and Appendix lA 3p . 

For open models (or LTB regions in which ([9]) allows for 
i?' > to hold everywhere) Ri can be any monotonously 
increasing function complying with Ri (0) — and R[> Q 
for all r. Evidently, the simplest choice in these cases is 



Ri = i?o r. 



(44) 



where i?o is an arbitrary constant characteristic length 
scale. The gauge ([i^ is a popular choice in the literature 
[l3 - [l7j , not only due to its simplicity, but because setting 
r = Ri/Ro allows us to regard radial dependence as a 
dependence on an fiducial value of an invariant quantity 
{R) that has a clear physical and geometric meaning. 
Also, the choice of Rq provides a physical length scale 
for the radial coordinate. 

However, other radial coordinate gauges are often used 
(fo r example setting M as radial coordinate, as in 0- 
[lOl|). And, as pointed before, ([44]) cannot be used in cases 
where a regular zero of R' occurs. For closed models, Ri 
must be selected so that it vanishes at both centers of 
symmetry and R'i has a common same order zero with 
J- = vl + E (to avoid the surface layer singularity asso- 
ciated with (|123|) . see [13]). Moreover, even if we use the 
parametrization associated with ([M]) . we are not forced 
to use the radial coordinate gauge for fixing Ri , it is still 
possible to keep this initial value function unspecified and 
use the gauge freedom to fix any one of the other initial 
value functions. 



VI. THE ANALYTIC SOLUTIONS IN TERMS 
OF AN INITIAL VALUE APPROACH. 



We rephrase in this section the analytic solutions ([5])- 
([7]) in terms of {iriqi, kqi}, which are related to the free 
parameters M and E in ([5])-([7]) by 



M = ruq, Rf 



E — kqi 



(45) 



A third free function (the "bang time", tbb) necessarily 
appears as an "integration constant" in both, the solu- 
tions of L in ([38]) and those of i? in ([3]). However, as 
we show below, if we rewrite ([5])~([7]) in the context of an 
initial value problem associated with niqi , kqi and L (con- 
sidering that Li — 1), the bang time i^b can be obtained 
as a function of mqi , kqi . 

Of course, other combinations of initial value functions, 
such as {niqi, Hqi} or {kqi, Hqi}, can also be considered, 
as ([38]) relates the three available functions. Since the 
fluctuations follow from the gradients of these functions 
from ()3ip . any two of {niqi, kqi, Hqi} forms an irreducible 
set of initial value functions. As we show in section XI, 
these initial value functions provide also the basic initial 
conditions for the evolution equations under a numeric 
approach. 



A. Parabolic models or regions: kqi — 

We express AI and R in terms of mqi and L from ([36] 
and ([45]) in ([5]), after re-arranging terms we get the fol- 
lowing closed analytic expression for L 



L 



1 + -y/2mqi C{t - ti) 



2/3 



(46) 



where we are only considering expanding configurations 
{L increases for t > ti). The bang time follows by setting 
L = and t = Ui, in ^ 



3^27 



Cti 



iT~Lqi 



(47) 



B. Hyperbolic models or regions: kqi < 

We obtain the following implicit solution of the form 
t — t{R) by eliminating the parameter rj from the equa- 
tion for i? in ([5]) and substituting in the equation for t: 



c{t~U^) ^ Zh{R), 



(48) 



where R — {E/M)R and Zh is the function 

u i-> Zh{u) ^ (2 + u)^/^ - arccosh(l + u). (49) 

We express then M, E and R in ([48]) in terms of ruqi , kqi 
and L from ^ and The resuh is 



Vi c{t - U) = Zh{xiL) -~ Zh{xi), 



(50) 
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where 



Vi 



(51) 



Setting t — and L = in (f50| and using (j49|) yields 
the bang time as a function of niqi and | fcgj | : 



Zh{xi) 



(52) 



C. Elliptic models or regions: kqi > 

In this case, the imphcit solution t — t{R) follows from 
eliminating from the first equation in ^ and substi- 
tuting in the second one. The resulting implicit solution 
has two branches, an "expanding" one (0 < 77 < tt with 
i? > 0) and a "collapsing" one (tt < 77 < 27r with i? < 0): 

\E\'^/'^ , ^ _ ( Ze{R) expanding phase 

^ ""'^[271- Ze{R) collapsing phase ' 



(53) 

where R = {\E\/M)R and is given by 

u ^ Zeiu) = arccos(l - m) - m^/^ (2 - u)^^^ . (54) 



Proceeding as in the hyperbolic case, we transform ((53 
into 



Ui c{t-ti)+ZeiXi) 



where 



Zf,[xiL) 



expanding phase 



2tt — Ze{xiL) collapsing phase 

(55) 



Vi = 



, 3/2 



''qi 



(56) 



It follows from ([55]) -([541) that L is restricted by < L < 
imax, where the maximal expansion is L^^^ = 2/xi = 
2mqi/kqi, characterized by i = and Hq — 0. 

Setting t = and L = in the expanding phase 
of (I55t yields a bang time function, while t = t^^n and 
i = in the collapsing phase yields the "crunch" time 
associated with the collapsing singularity. The maximal 
expansion time follows by substituting t = i,„-,x and L — 
in either branch of (j55p . These times are given by 



C^bb = ct 



Ze{xi) 

Vi 



(57a) 



TT TT — Zp(Xi) , , , 

= cibb + - = cU + (57b) 

27r 27r — Ze(Xi) ^ 

ct^^u = ctbb H =cti-\ . (57c) 

Vi Vi 

Notice that (in general) t^^^ = t^^^^{r) and t<,„„ = tcoii('^): 
so neither one coincides with a ^T[t\ hypersurface (like 



^bb(''))- For every comoving observer r = const., the time 
evolution is contained in the range tbb(?') < t < icoii('')- 

Rewriting the analytic solutions ©-([T]) as in and 
as the implicit forms ([50]) and ([55]) is very useful, as it al- 
lows us to relate (by implicit radial derivation) the gradi- 
ents of L (and thus of R) with gradients of the intial value 
functions rriiq, kiq (see section 11). Notice that the time 
splitting associated with m° and the scaling of R with Ri 
in L imposes on LTB models a constraint between tbb 



and niqi, kqi, since (|47p and (|52p imply that any choice 
of these functions uniquely determines the "age" ti — tbb 
of the initial slice t = ti for any comoving observer. Also, 
the initial slice t = ti in elliptic models always occurs in 
the expanding phase [2^ . 



VII. CURVATURE SINGULARITIES. 

The new variables allow us to express the coordinate 
locus of curvature singularities exclusively in terms of L 
and r by means of (gl]), (gS]) and (gS]). Notice that for 
reasonable initial value functions (bounded and contin- 



iq, kq, Hq diverge as L 



0, 



uous), all scalars Aq 
whereas local scalars A = m, k, H can also diverge if 
r — (even if L > 0). Following the standard criteria 
that define a curvature singularity and considering (jA4l) , 
we can identify two known possible singular surfaces 

L{t,r) — 0, central singularity (58a) 
T{t, r) — shell crossing singularity. (58b) 

Notice that if F > 0, then all scalars A and Aq only 
diverge at the central singularity L = 0, but if F — >■ 
0, then all the relative fluctuations diverge (with 
Aq 7^ 0), SO that local scalars A diverge while their quasi- 
local duals Aq remain bounded. This is an obviously 
unphysical effect of shell crossings that must be avoided. 

It is important to remark that characterizing the cur- 
vature singularities as in (|58ap and (|58bp represents an 
improvement in clarity over the conventional variables. 
In these latter variables "i? = 0" can be either a regular 
symmetry center, or a curvature singularity. Likewise, 
"i?' = 0" can be either a shell crossing singularity or a 
regular "turning value" r = rtv associated with a closed 
model. There is no such ambiguity in our parametriza- 
tion, since L > holds at a symmetry center and F > 
holds at the turning value. 

In order to test the avoidance of (j58bp we need to 
compute F. For this purpose, we derive both sides of 
the solutions ([50]) and (|55|) . then we use (I5T|) and 

L'/L — {\~T)R[/ Ri to ehminate m^j, A:^ and L' in terms 

of 5^^\ 5^^ and F. The result is 



Parabolic models or regions. 



F = 1 + - 



L3/2- 



(59) 
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Hyperbolic and elliptic models or regions. 



A. Parabolic models/regions. 



r = 1 + 3(<5^ 



(m) 



<5f )) ( 1 



T~Lq 
Ho 



c(fe) 



(60) 



where Hq and Hqi follow from (|38| . while c{t — ti) 
is given by and ([55]). 

Inserting (|59p or (|60l) (depending on the kinematic class) 
into (|4T|) . (l42t and (|43)) leads to closed analytic expres- 
sions for all these scaling laws in terms of L and ini- 
tial value functions in the context of the initial value 
approach described in previous sections. The functional 
forms of r and T-Lq in ((38|) and ((59)) allow us to obtain 
analytic expressions for all the scaling laws of covariant 
scalars (as functions of L and initial value functions) . 

From dH), dSI]), (go]) and (l42al) applied to A = m, con- 
dition dH) for standard regularity that avoids (|58b|) now 
reads 

r > V (c<, r) such that L > (61) 

The fulfillment of this condition is examined in the fol- 
lowing section. 



VIII. THE HELLABY-LAKE CONDITIONS. 

We obtain in this section the Hellaby-Lake conditions 
JTUll-dlll) that fulfill the regularity condition dH]) [2i, .25, 
|39| . but now as restrictions on the initial value functions. 
These conditions guarantee absence of shell crossings for 
all t. Notice from © and (PH)) that avoidance of this 
singularity already excludes a surface layer singularity 
associated with (|123l) . 

In order to examine the Hellaby-Lake conditions in the 
initial value parametrization, we need to relate the gra- 
dients M' and E' appearing in pIH) - ([Tl]) with the initial 
fluctuations by specializing to t = ti 



M' 



3^ 



1 



El. 
'e 



(62) 



We will also need an expression for the radial gradient t'^^^ 
in terms of our initial value functions. This expression 
follows by differentiating both sides of (l47l) and ((52)) with 
respect to r: 



3i?^/ Ri 



2 6, 



(m) 



3^2mq 



ct' 



?>R[/ Ri 



\ i 

T~Lai 



c(ti - ibb) 6, 



-Am) 



parabolic 
(63a) 

2 * 

hyperbolic and elliptic (63b) 

where in ()63bp is given by (|52]) with Zh{xi) or Ze{xi), 
respectively, for hyperbolic and elliptic models. As men- 
tioned before, Ri can be always be prescribed as a radial 
coordinate gauge. 



From (HD), if L ^ then F -S'f'^/L^/'^, while F -> 

1 + ij'™' for L — > oo. Considering ()42ap . the necessary 
and sufficient condition for ()6ip to hold is simply 



1 < -5^™^ < 0, 



(64) 



which from ()3ip implies m'^^ < 0. It is evident from ()62p 
and ()63al) that dMj) implies < and M' > 0, and 

so it is equivalent to the Hellaby-Lake conditions ([TU)). 



B. Hyperbolic models/regions. 

For a fixed r all initial value functions are finite, con- 
sider then ti > ti > tbb but ti « ttt, so that L « and 



"Hg — > oo, and thus (IMl) becomes 



SHq 
T~Lai 



— ^SHa 



Am) 



ctL 



Ak) 



Hqi C{ti -tl) [6. 



Am) 



2 ' 



(65) 



where we used (j63bP . Now consider t — oo so that L ^ 
oo for all r. From (l38l) and (l50l) we have in this limit 



"Hg — > and 



HqC{t-t,) « 1 + 



Zh.(xj) -|-ln(2xjL) - 2 



hence F in (l60l) becomes 



1 



|<5f' +0(L-2). 



-0(L-2), (66) 



(67) 



Bearing in mind, from ^ and (|A6p . that i?'; > for 
hyperbolic models/regions and that ([CT)) must hold as 
^1 — ^ ibb in (|55)) . and since the scaling law (|42al) and 
(|65 )) - (|67)) must hold for all t and r, the conditions for 
m are [13 



ci' . < 0, 



Ak) 



> 



<5f"') > -1. 



(68) 



where t'^^^ is given in terms of our initial value functions 
by (I63bp . Notice, from (I42bp . that all regular hyperbolic 
models comply with S'^^^ > —2/3 for all their time evolu- 
tion. From (|34| and (|62)) . it is evident that (|68)) are com- 
pletely equivalent to the Hellaby-Lake conditions ()10l) . 
Since these conditions are neces sary and sufficient, (|68)) 
are also necessary and sufficient [Ij, IH, [s^ . 



C. Elliptic models/regions. 

We look again at F in (jSO)) . but now considering that 
T-Lq changes sign (from positive to negative) at some t = 



10 



^max where L = L,nax = '^rriqi/kqi = 2/xi, so that < L < 
Linax and we have the presence of a second (coUapsing) 
singularity as L ^ when t — ^ t^^„, with t,„-,x and ^con 
given by (j57p . Near the maximal expansion L — ;> i^ax 



and so Tiq — >■ 0, while c(t — ti 



c(t„ 



ti). Thus, we 



have in this limit the following necessary (not sufficient) 
condition for (1611): 



r = l + 3((5|™^-(5fO >o 



(69) 



Near the initial singularity (bang) we have "Hq — > oo as 
t thb, while near the collapsing singularity we have 
"Hg— >— ooast— )■ tcoii with t^^u given by ([57|) . We 
consider t — ti such that U > ti > but ti ttt, as 
well as t = t2 such that tcon > t2 > U with ^2 ~ ^cou- For 
these times F takes the form 



which, because of (l42c)) . implies that 5^™) - (3/2)(5('=) > 
necessarily holds for all times if (|6T1) holds. Since condi- 
tion ([74|) implies (|69|) . then the latter (also necessary but 
not sufficient) follows from (|72l) too. 

It is straightforward to show from ([3T|) . ((57l) . (|62|) . 
(|63bp and ([73l) that conditions (|72|) are equivalent to 
the Hellaby-Lake conditions (fTTj). We have expressed 

these conditions in terms of rriqi, kqi, i5j^"'\ but it 
is straightforward to rewrite them in terms of other sets 
initial value functions, like m^i, T-Lqi, 5'f^\ S^^K Because 



of (jA6P and (j40| , we do not need to express (|72p in terms 
of the ± signs for R' , as it is understood that the signs 
of t'^^^ and i'^ij must be, respectively, the same and the 
opposite of the sign of i?- (which because of (|A6P will be 
the same as the sign of R'). 



SHq 
T~Lai 



ct' 



(fe) 



3 



.(fc) 



iR'jR^ 



+ c(ti - ibb) 5, 



(70) 



« 3^ 

T~Lqi 



Hq, c{t2 - u) [ s; 

cti 



(m) 



2 ' 



(71) 



Since ([70]) must hold as ti — > i^b and ([7T|) must hold as 
^2 ~^ ^coiu each expression furnishes a second and third 
necessary condition for (|5T|). Since both conditions that 
follow from (ITO)) and (|7T]) must hold for all t in the range 
tbb{r) < t < t^au{r) for all r, these two conditions taken 
together are then the necessary and sufficient conditions 
for the fulfillment of (j6ip . which can be written as 



ctL 



m'jR., 

with 
ct' 



< 0, 



ct' 



iR'jR, 



>0, 



> -1, 



3i?^/ Ri 



girn) _ £^(fc) 



c(icolI ^bb) + 



3-R^/i?, 



(72) 



(73) 



where we used c(tooii — ti) — c(tooii — ^bb ) — c(i i — tbb ) , with 
ibb and t^oi, given by (|57l) . and the equal sign holds only 
at a symmetry center. The fact that one of the three con- 
ditions in (17^ to avoid shell crossings in elliptic models 
is basically a sign condition on the gradient t'^^^^ has not 
been, apparently, noticed in the extensive literature that 
uses the conventional variables. The reader is requested 
to compare the simplicity and elegance of (|72p with the 
cumbersome form (jlip . 

Conditions ()72p imply the following necessary (but not 
sufficient) condition for ([BT]) 



>o. 



(74) 



IX. OMEGA AND HUBBLE EXPANSION 
PARAMETERS. 

The fact that Hq behaves in ([55)1 as a Hubble scalar of a 
FLRW dust universe, has motivated some authors [ll[ll| 
dealing with cosmological and astrophysical applications 
to define an "Omega" parameter analogous to that of a 
FLRW dust cosmology 



KPq_ 

mi 



2mq 
Ti? 



2m„ 



2ma 



2ma 



so that 



n - 1 



2mo 



kqiL 



kqiL 



2'mq kq 27Tlqi 



kqiL 



(75) 



(76) 



where we use the symbol fJ, and not fi^, because this 
quantity is not the quasi-local dual of the scalar 2m/H^ 
under the definition (|24)l (see section XIII) . The following 
FLRW-like scaling laws hold for Ct and 17 — 1 : 



Qi 



- {n, -i)L 



(n, - i)L 



- {n, -i)L 



(77) 

while the Hamiltonian constraint psp can be rewritten 
as an analogue of the FLRW Hubble expansion 



(78) 



The relative fluctuation P5)) of the quasi-local Hubble 
factor, takes the appealing form 



(79) 



which can be given in terms of metric functions L, F 
and initial value functions Cli, Jj-"*-* and 5^^^ by means of 
g^ . (I42b)) . O and 

We can identify di in ([77)) . (1751) and ([7^ as analogous 
to an r— dependent FLRW Omega factor at a fiducial 
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cosmic time ti 



stitute 



2m 



qi 



2mn 



1/2 



2m„ 



(80) 

Notice that the sign of f2i — 1 in any region of r deter- 
mines the sign of 57 — 1 for aU times, which is the same 
as with kqi and kq. Hence, Cl — 1 behaves as the FLRW 
Omega factor for spatial curvature, so that 0.^ — 1 deter- 
mines the kinematic class: parabolic if fi^ = 1, elliptic if 
r^i > 1 and hyperbolic if < f2i < 1. Also, irrespective 
of the kinematic class we have for every LTB model or 
region: — > 1 as L — > (near the central curvature 
singularity). For all hyperbolic models or regions VL is 
bounded between and 1 for all choices of 17,;, T-Lqi (or 
Wgi, kqi), hence: £7 — )■ as i — > oo, whereas for elliptic 



models or regions il oo as L L^^^ because T-Lq 



0. 



Too 



Til - 1) 



(81) 



in (|50)) and ([55]). Since tli, Hqi are roughly equivalent 
to inhomogeneous generalization of fiducial Omega and 
Hubble factors, parametrizing the analytic solutions with 
these initial value functions could be more intuitive than 
doing it with the density and spatial curvature profiles 
niqi, kqi. The solutions (150]) and ([55]) become 

• Hyperbolic models or regions: ilj — 1 < 0. 



c{t - U) 



n 



qi 



• Elliptic models or regions: f2j — 1 > 0. 



(82) 



A. Analytic solutions in terms of Cli 

It is straightforward to parametrize the analytic solu- 
tions in terms of fi^, Tiqi- For the parabolic case we have 



W — Wi expanding phase 
■Hq^ c{t -U) = { ^^^^^^ _ ^^_3/2 _ ^ _ (83) 
collapsing phase 



Cti — 1, and thus we just obtain with yj2mqi = Tiqi- 

For the hyperbolic and elliptic models or regions we sub- where the functions W and Wi are 



Hyperbolic models or regions 

W 
W, 



Elliptic models or regions 



+ (1 - h,)L 



1/2 



f2L 

— arccosh — — h 1 — 2L 

2(1-17^)3/2 



2(i-Oi)3/2 



-arccosh 



a-l 



W = — arccos — + 1 - 2L 

2(17,-1)3/2 

w ^» /2 \ 1 
Wi — — arccos { 1 — . 

2(r!,- 1)3/2 y^i^ J 



1/2 



Li/2 



r 



(84b) 



(85a) 
(85b) 



Setting t = and L = in (j82p and in the expanding 
phase of ([83]) yields the bang time for hyperbolic and 
elliptic configurations 



while the maximal expansion and collapse times, t — t„ 
and t — i„„„, given by (1571) are 



^max = ctb 



^coll 



7ri7,; 



2nq^[n^-l]y^' 

nCli 

^H,,[f],-l]3/2" 



(87a) 
(87b) 



Cti 



W^ 



(86) 



We can also parametrize F, as well as the Hellaby-Lake 
conditions, in terms of i7i, Tiqi and their gradients. For 
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the parabolic case, we substitute into ((59)) and (|64] 



a constraint that will not be satisfied unless mqi and mqj 
have very restricted forms. For a hyperbolic model we 
rewrite (l50l) as 



while for hyperbolic and elliptic models or regions we 
substitute into ((501), dSH]) and ^ equation ^ plus 



F{\kqA) 



Zh{.a\kqi\) - Zh{l3\kqi\) 



qi 



13/2 



c{t, -U), (93) 



R^_ 

Ri 



9i 



T~Lqi 



where Zh is given by (j49|) and 



9* 



(89a) 
(89b) 



1 



2/3 1/3' 



1 

m, 



(94) 



Since a unique LTB model is determined by prescribing 



niqi and \kqi\, and we assume to. 



9« 



and 



known, then 



where i?, can be specified as a choice of radial coordinate. dSSl) becomes a constraint to find the missing initial value 

function \kqi\ that determines the LTB model. Solving 
this constraint for known to^^, mqj and c{tj —U) can only 
be done numerically, but by looking at the properties of 
the function F{\kqi\) we can find the conditions for the 
existence of this solution, as it was done in 0, • 

Considering that in an expanding model we have mqi > 
mqj for tj > ti, we have a > /3 in general, and since Zh is 



The initial value functions 17,, Tiqi can also be used 
in the analytic solutions given in the forms ©-([T]), in 
terms of the parameter 77. For this purpose, we eliminate 
AI and E by means of 

^ _ 2Af 2 _ 2M + ER, 

where Ri can be prescribed as a radial coordinate gauge. 
The bang time is then given by 



monotonously increasing, then the left hand side of ([93 
is positive (like the right hand side). We also have 



F 



X. MAPPING OF DENSITY TO DENSITY OR 
VELOCITY TO VELOCITY PROFILES. 

Krasihski and Hellaby 0, have examined an inter- 
esting proposal: the possibility that arbitrary radial pro- 
files of density (or velocity or density and velocity) given 
at two different cosmic times for a finite radial range can 
be fit (or "mapped" ) by a unique LTB model. The initial 
value parametrization of the analytic solutions is ideally 
suited to examine this issue, which can also be under- 
stood in terms of the equivalence of boundary and initial 
conditions to solve the evolution equations of sections 
XI and XIII. Our aim here is to illustrate how the ini- 
tial value parametrization can be helpful, hence we only 
examine in detail the mapping between the quasi-local 
density profiles by a hyperbolic model, since all other 
similar mappings can be handled in the same manner. 

Consider two radial profiles of m^. Since t — ti is 
arbitrary, then mqi can be one of the profiles, while we 



{\kqA)^^{c.'''-P''')+0{\kq^l 

a- 13 



qi\ 



F{\kq^) 



/3|fc, 



0{\kqA~"^). 



qi I 



- 0, 
(95a) 
00, 
(95b) 



which, together with similar expansions of dF/d\kqi\^ 
show that _F(|fcgi|) is monotonously decreasing and 
bounded by 



-1/2 -1/2 



Q<F{\kq,\)<^(a^'^~P''^)=^ 

(96) 

Hence, considering (|9T1) . a unique solution of (|93l) must 
exist if 



Z \ mqj J 



1/2 



= L 



■i/2 



(97) 



take 



TOq {ctj , r) for some tj > ti as the second 



profile. If mqi and mqj are known, then 



Lj{r) — L{ctj, r) 



1/3 



(91) 



follows from ([BT)) (remember that Li — 1). For a 
parabolic model, we have from (j46p 



c(tj -t,) = 



3^27 



1/2 



(92) 



which is the same result found in Comparison with 
(j46| shows that the density-density mapping requires the 
LTB model that maps mqi to mqj to expand faster than 
a parabolic model. For the elliptic case we follow the 
same method: using the solution (|55|) (in the expanding 
and collapsing phases) as a constraint to find the missing 
initial value function kqi . If the local density profiles to^ 
and TOj are known (instead of mqi and m„i), then we 

ti r 

{Ri Lj 



have to use the scaling law ()41a|) with F given by (|40] 
3to 



^2 



iLj Tj RiR'i 



'en 

3 r3\ / 



(98) 
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which yields Lj by integration (as and ruj are known 
and Ri can be fixed as a coordinate gauge) . It is straight- 
forward to show that the analytic solutions for the hy- 
perbolic class yield a unique LTB model mapping these 
profiles if (HH) holds. 

Krasinski and Hellaby also examined the mapping of 
profiles of velocities defined as w = R/Ri — LHq. To deal 
with this case we eliminate kg so 
that the initial value functions are m^j and T-Lqi- Then, 
from the Hamiltonian constraint (1381) we obtain 



2mq,Rf 



2mq,Rl '' 



(99) 



which inserted in the solutions ((50)) and ([55]) yields the 
constraints needed to obtain the missing initial value 
function rUqi. We can also obtain the LTB model that 
results from assuming that the profiles fli and flj are 
known. In this case we obtain Lj from (j77|) 



n, 1 - rij 



(100) 



which substituted in the solutions and ([55)) yields 
the missing initial value function Jiqi by means of the 
constraint 



XI. EVOLUTION EQUATIONS FOR A 
NUMERICAL TREATMENT. 



The analytic solutions, either in terms of the con- 
ventional variables M, E, cttb, R or the initial value 
parametrization rriqi, kqi, Ri, L that we have presented, 
are mostly useful for qualitative work (as illustrated in 
the previous section). We remark that these solutions 
are either parametric or implicit, and thus a pure ana- 
lytic framework based on them is necessary limited. For 
concrete applications in models of less idealized cosmo- 
logical inhomogeneities and observations there is no other 
way but to consider a numeric framework. 

The 1-1-3 fluid flow system ([T5)) ~(P(7 1) discussed in sec- 
tion IV provides a fully covariant description and can be 
used to determine the dynamics of LTB models by means 
of numeric techniques. However, the associated spacelike 
constraints (I19|) are partial differential equation on r that 
are not easy to solve, making it difficult (in general) to 
decouple the time and radial derivatives. As we show in 
this section, the quasi-local scalars and their fiuctuations 
(in the initial value parametrization) provide more con- 
venient variables for a numeric treatment of LTB models. 



A. Quasi— local evolution equations. 



"Hqi — 



c{tj -u) 



(101) 



where W and Wi take the forms (|84al) - (j84bp or (I85ap - 
(|85b[) with L given by (|100|) for hyperbolic or elliptic 
models or regions. Evidently, obtaining the missing ini- 
tial value function by means of constraints like (j93p and 
(jlOll) does not guarantee the fulfillment of regularity con- 
ditions (like absence of shell crossings). This has to be 
verified independently by means of the Hellaby-Lake con- 
ditions once the LTB model mapping the profiles has 
been found. 

The issue of mapping two profiles can also be under- 
stood in terms of the systems of evolution equations that 
we derive in the following section. These systems are 
partial differential equations, hence initial and boundary 
conditions are given as functions of r (or radial profiles). 
The selection of two initial value functions is sufficient 
to determine a unique solution. These functions can be, 
either niqi, T-Lqi or mqi,kqi or Vli, T-Lqi. Hence if we know 
only one of these functions (evaluated ai t — ti) and the 
same function at t = tj , then we may use one of the avail- 
able constraints to find the second initial value function, 
and thus guarantee (if the constraint has a solution) that 
the evolution system yields a unique LTB model. We 
have proceeded in this section (following 7, 10]) by us- 
ing the analytic solutions as the constraints that link the 
boundary and initial conditions through L. The results 
are the same as those of 0, [13, but the initial value 
parametrization of the solutions makes the whole proce- 
dure a lot more natural and straightforward. 



An alternative system to the evolution equations ([18)) - 
([^01) follows by considering a 1-1-3 ffuid fiow evolu- 
tion equations for the quasi-local scalars (|35p . involv- 
ing mg, Tiq and their corresponding relative fluctuations 
^(m)^ This system can be obtained by eliminating 

p, O, S and £ in ([181) in terms of the new variables by 
means of ([5^ and using ([5T)) to deal with the constraints 
(dD). Considering the notation H = 6/3, to = 

Kp/3 and k = ^TZ/Q, this framework leads to the follow- 
ing evolution equations 



-SmqTiq, 



Tiq — —Tiq — niq, 



-3(l + (5("))Hg(5(^), 



(102a) 
(102b) 
(102c) 



-(l + <5(«))H,<5(«) + ^(<5(«)-<5(™)), 

Tiq 



(102d) 

while the spacelike constrains (IT^ become the equations 
defining and S^'^^ through ([27a)) and ((311), and the 
Hamiltonian constraint ([^0)) reduces to ([381) . This is the 
particular case of pure dust in [so', 'sT] . 

The system ([1021) is equivalent to ([T5)) . but it is much 
more practical and easier to use in a numerical treat- 
ment, as it is not necessary to solve any radial differen- 
tial equation as a precondition to solve it. In practice, 
we can handle it as a system of non-linear autonomous 
ordinary differential equations, where r enters as a pa- 
rameter (see Appendix B of [231). Notice that once we 
have solved (I102p all local scalars in ([T7)) can be deter- 
mined by means of ((23) and (H51) . 
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diverges as Hq — > 0, as the maximal expan- 
sion ct = ct^.^^ is reached in elhptic models or regions 
fsee I A 4p . it is better to use in these cases the alterna- 
tive variable S = -HqS'-'^') (see This leads to the 
alternative system 



Hq — 



sir' 



-imq'Hq, 

-Hi - mq, 
-3(1 + ^(") 



E = -2Y.nq 



ruq S'^" 



(103a) 
(103b) 
(103c) 
(103d) 



which is entirely equivalent to ()102p . and has the advan- 
tage (besides its use for recollapsing configurations) that 
S is the scalar associated with the shear tensor by ([T5|) . 

We may supplement either (|102[) or (|103p with the fol- 
lowing extra pair of differential equations 



L 



L Ha 



t = iV'HqS^'^'^ = -3rs, 



(104a) 
(104b) 



which follow directly from ([25]) and (gO]). By includ- 
ing these extra equations we can use the systems above 
to determine numerically the metric functions in (1391) . 
In fact, if we substitute the scaling laws ([37)) -(l38 |) . (|411) . 
(|42| and (|43|) into ()104p we can work only with these two 
evolution equations given as 



L 

r 



[2mgi — kqiL\ 
3 2m,, [<5,|") 



1/2 



(105a) 



i-r]-fc,,£[jf^ + |(i-r)] 

[2m,,-fc,,i]i/2 



(105b) 



whose solution {L and F) allows us to compute all scalars 
by means of these scaling laws. However, from a numeric 
point of view it might be more convenient and practical 
to work with or pU5|) supplemented by pMl) than 

with (fT05)) . 

Solving any one of the systems of evolution equations 
described above is needed to find radial null geodesies 
and their corresponding red shift factor, z, which are in 
turn needed for computing observable quantities (lumi- 
nosity distance, etc). Given a numerical solution of either 
system, radial geodesies ct = cti^(r) in the past null cone 
and z follow from solving 



dr 



cdtAf(r) 
dr 



ln(l + z) 



N 



(1 + ^)1/2 



Ln Tjv R'i 

[i-kq,E:iY/^' 

(106a) 



N 



(1+^)1/2 
Ljq Fat HqN (1 



(106b) 



where the subindex means evaluation at the null curve 
ct = ctAr(r), i.e. = A{ctN{r),r). Since the systems 



(|102p or (|103|) or (|105l) can be handled in practice as 
constrained ordinary differential equations, solving nu- 
merically (|106|) can be achieved with simple numerical 
techniques. 



B. Initial conditions. 



An obvious choice of initial value functions for the sys- 
tems pMt - (fT04)) or ([T03l) -(IT04 | is 



n 



qi I 



j(m) _ ^T-'qi/'l^qi 

* ZR'JRi 



SR[/ Ri 



iR[/ Ri 



(107) 



where we used (|3ip and Ri can be prescribed as a radial 
coordinate gauge. 

Since kqi determines the solutions of (p8|) (or ^ for 
R = RiL), this initial value function is closely related to 
the kinematic evolution of the models (or specific regions 
of them). Hence, it is more practical and intuitive to 
choose instead of (|107l) the initial value functions 



^qi , 



e(m) 



qi 



SR'jRr 



k'q^/k. 



qi 



' 3R[/Ri 
(108) 

These initial conditions are sufficient to solve ()105p . but 
to solve pII^ - ([TM)) or ^U^-i^U^ we also need "H,, and 



or Si, which readily follow from (|38|) and (j43 



Ho 



[2mq,-kq,]'/^, 



or 



2[2mqi - kqi] 



2[2m,,-A:,,]i/2 ' 



(109) 



The main advantage of using (llOSp as initial conditions 
is the fact that they are exactly the same initial value 
functions that we employed in the parametrization of the 
analytic solutions and in the forms of F in (l5n|) - (|M| . and 
the Hellaby-Lake conditions. Hence, given a choice of 
these functions, a time evolution free from shell crossings 
can be immediately tested by means of (|64p . ([SS)) and 
dZl]). 

The evolution equations (|103l) and have already been 
used for a dynamical systems approach to LTB mod- 
els [29], while a suitable generalization of (|102p was 
employed for a numeric study of cosmological models 
endowed with the LTB metric but with an energy- 
momentum tensor of an anisotropic fluid [31]. 

A comprehensive numeric study of LTB models by 
means of the systems (|102p . (|103p or (llOSp is beyond 
the scope of this article, and thus is presently under ex- 
amination in a separate work. Our purpose has been 
to show how these systems, constructed with quasi-local 
variables, are potentially promising for this purpose. 
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XII. GAUGE INVARIANT AND COVARIANT 
SPHERICAL PERTURBATIONS ON A FLRW 
BACKGROUND. 



The connection between the system (|102p and a per- 
turbation formaHsm on a FLRW background is clearly 
suggested by the fact that the evolution equations (|102ap 
and (|102b|) for the quasi-local scalars niq, Hq are for- 
mally identical to the energy balance and Raychaudhuri 
equations of a dust FLRW cosmology. Also, from 
(IMt. (I37l) and (ESI, we have L 



—2kq'Hq, which is the 
same evolution equation for '^TZ in a FLRW spacetime, 
while the 5*^"^^ defined by ([?T|) and its relation with the 
scalars in ([5^ clearly highlights a sort of perturbation 
definition. We provide in this section a rigorous charac- 
terization of this intuitive resemblance. 

A perturbation formalism between an idealized space- 
time, S (a dust FLRW cosmology), and an "lumpy" uni- 
verse model, S (a dust LTB model), follows by defining 
a "background model" in S, constructed by objects (in 
S) that are the images of a suitable map $ between ob- 
jects in S to objects in S [s^. Though, perturbations 
in general are not uniquely defined by such an abstract 
map because of the ambiguity of the gauge freedom 
associated with choices of coordinates and hypersurfaces. 
This fact requires carefully defining perturbed variables 
that are "gauge-invariant" However, in the specific 
case that concerns us we face a simpler task, as both the 
idealized FLRW and the perturbed LTB spacetimes are 
LRS (locally rotationally symmetric) and can be com- 
pletely described by scalar functions [s^. Hence, $ can 
be a map of scalars to scalars, while the fact that both 
spacetimes are spherically symmetric and are given in the 
same normal geodesic coordinate (or frame) representa- 
tion, greatly suppresses most of the gauge freedom that 
we would find for a general 5'.(48j 

The map $ can be defined rigorously as follows. Let 
X and X be, respectively, the sets of smooth integrable 
scalar functions in S and S, then for all covariant FLRW 
scalars A G X (we denote FLRW objects with an over- 
bar) the map 



^ -.X ^ X, $(A) ^ AqeX, 



(110) 



defines a "background model" (associated to a FLRW 
cosmology) in LTB models through the quasi-local 
scalars Aq (which are LTB objects satisfying FLRW dy- 
namics). Given their common normal geodesic represen- 
tation and considering the perturbations 



A-$(A) 
HA) 



(111) 



associated with ()110|) . LTB models characterized by the 
metric ([T]) and source ([2]) can be considered then as spher- 
ical non-linear "perturbed" dust FLRW cosmologies in 
the "comoving" gauge. 

The perturbation scheme that we described fits nat- 
urally to the gauge invariant and covariant (GIG) ap- 
proach of Dunsbury, Ellis and Bruni [s^ [s^ . Following 



these authors, a perturbation scheme on FLRW cosmolo- 
gies is covariant if the "lumpy" model S is described 
by variables defined in the framework of the 1+3 fluid 
flow variables associated with the system ^TE\\ . Although 
our description of LTB spacetimes is not based on these 
scalars (the representation (IT71) ') but on the quasi-local 
representation (1351) . it is still a covariant description be- 
cause rriq and Hq , are covariant scalars by virtue of their 
connection with the invariants M, R and their deriva- 
tives in (|28l) . Hence, the formalism based on niq and Hq 
would also be covariant. 

As commented by Ellis and Bruni [s^l , by virtue of the 
Stewart- Walker gauge invariance lemma [35[, all covari- 
ant objects in S that would vanish in the background S (a 
FLRW cosmology in this case) are gauge invariant (GI), 
to all orders, and also in the usual sense (as in [33]). The 
tensorial quantities in LTB models that vanish for a dust 
FLRW cosmology in the l-|-3 formalism are the shear and 
electric Weyl tensors tr"'' and given by (IT51) and 
in terms of the scalar functions E and £, which from ([501 . 



can be written in terms of the fluctuations m 



and 



H — Hq. From (PT|) . it is evident that these fluctuations 
and S'^'^\ as well as the radial gradients m'q and "H^ are 
all "flrst order" quantities (in the perturbation scheme) 
that are GI to all orders, though niq, Hq are not (which 
is expected because these are "zero order" background 
variables). As a consequence, the fluid flow dynamics of 
LTB models in the quasi-local scalar representation ([55)1 
also provides a rigorous characterization of LTB models 
as spherical, non-linear GIG perturbations on a FLRW 
background. 

We note that spatial curvature k — and its quasi- 

local dual kq — ^TZq/Q do not appear in the dynamic 
equations p02ap - (ll02dl) . though they can be used when 
specifying initial conditions as in (llOSp (see also the ap- 
pendices of [111). Spatial curvature is GI only if the 
FLRW cosmology S is spatially flat, but its associated 
variables 5^^'^ and k'q are GI. In fact, from ([55]) . (|571) and 
([55)1 we can always eliminate either one of m,, T-Lq or 
g(m) ^ gCH) terms of kq and S^'^\ and construct a sys- 
tem of evolution equations equivalent to (I102ap - (ll02dl) . 
but describing the dynamics in terms of spatial curvature 
kq and its perturbation S^''\ 

Given the correspondence between p02p and the dy- 
namics of non-linear perturbations on a FLRW back- 
ground, it is important to examine the connection with 
linear perturbation theory of dust sources. For this pur- 
pose we derive a second order equation for the density 
perturbation J^"*) by differentiating both sides of (|102ap 
and use the remaining equations (jl02bp - (|102dp to elimi- 
nate all derivatives except 5'--^^ and S^"^K We obtain 



sir' 



^^+^nqSi-^^lpqSi-^ (H-^(")) 



= 0, 
(112) 



which is an exact non-linear equation for (a sim- 

ilar equation was obtained in |41j]). For near homoge- 
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neous conditions, as assumed in linear perturbations with 
"small" perturbations \S'^"^'> \ < 1, (|112p reduces to 



2nq (5(' 



Pq 



sir 



0, 



(113) 



an equation that is formally identical to the evolution 
equation for linear density perturbations of a dust source 
around a FLRW background (characterized by pq, T-Lq) in 
the comoving gauge [38| , which for dust is a synchronous 
gauge as well. 



XIII. A THEORETICAL CONTEXT FOR THE 
OMEGA PARAMETER. 



parameters of a FRLW dust spacetime. As a consequence 
Cl is not a background variable in the perturbation for- 
malism, though its gradient defined as 

A = = (5(™)-2(5(^) = {l-d) (Jf")-,^^'^)), (116) 

oR / R 

is a GIC perturbation in this formalism (since 5^™^ and 
(5^^-' are), though we have A ^ [Vl — Q)/VL in general 
(this relation only holds in the linear limit). On the other 
hand, it is hard to find a theoretical context for Hcs, 
besides their resemblance to FLRW quantities and their 
utility in computations. 



The parameter (l that we discussed in section IX has 
been introduced in the literature [H, [l^ (together with 
T-Lq) as useful ansatzes justified by their resemblance to 
the corresponding FLRW Hubble and Omega parame- 
ters. Another generalization of the FLRW Hubble factor, 
suggested by Moffat and Tartarsky [i^ and used in var- 
ious articles follows by defining two expansion fac- 
tors, a radial and an azimuthal one, which can be given 
in terms of Hq and J^^-* as: 



Ha 
Hi 



R 
R 



L 
L 



Ha 



(114a) 



L 



(114b) 



Evidently, and -ff|| correspond, respectively, to ex- 
pansion factors of proper lengths in the direction orthog- 
onal and parallel to radial rays orthogonal to the orbits 
of S0(3). Moffat and Tartarsky define from these expan- 
sion factors an "effective" Hubble parameter as 



Hi 



Hi 



l~Lql~L 



(115) 



so that an Omega parameter follows as fieff — up/ {iH^^). 
However, both pairs "Hq, and H^s, fieS are basically 
useful quantities in the qualitative or numeric application 
of LTB models, since the proper local covariant general- 
ization of the FLRW Hubble parameter to LTB models 
is neither "Hq, nor ileff, but the expression given by equa- 
tion (41) of [1^. Nevertheless, it is still interesting to 
discuss the theoretical assumptions underlying these ex- 
pressions. 

While both T-Lq and TJeff are covariant quantities, the 
theoretical context for T-Lq may be easier to justify, as this 
scalar is the quasi-local dual of the fiuid flow expansion 
scalar H = 6/3, and thus it is a GIC background variable 
in the perturbation formalism discussed in section XII (it 
is the image of the FLRW Hubble parameter under the 
map (IllOp ). However, Ti,^ is not the quasi-local dual of 
"H , since the deflnition ^ implies that {U'^)q ^ [UqY, 
hence Ti.^ (and thus) Cl are not images under the per- 
turbation map (I110[) of the Omega and squared Hubble 



A. Evolution equations in terms of Cl and A. 



Since and T-Lg in (1751) and (iTSll are constructed from 



TTiq and kq , both quantities can easily be computed from 
a numeric solution of the systems (|102l) or (I103p . The 
initial value functions niqi and kqi in (I107p . (|108p and 
(I109P can be given in terms of Cli and Hqi from ((5T|) . 



while sl™'' and S}"^' follow from (l89 



(fc) 



However, given the practical utility of and the con- 
nection of its gradient A to the perturbation formalism 
of section XII, it is still useful to construct a system of 
evolution equations that include evolution laws for these 
quantities. Rewriting the system (|102p in terms of and 
A yields the following dimensionless system 



;F = -HM 1 



dA 
dr 



2 



H 



S^^'> + A ) O - ( 1 3A 4<5(^) 



(117a) 
(117b) 

(117c) 



= -H 



1 + <5(«)U(«) + ^(<5(«)+aU 



(117d) 



where we have introduced the dimensionless variables 



H; 



Hq 



HQc{t - U), 



(118) 



with Hq an inverse length scale (cm^^), which can be 
identified as the cosmological Hubble scale at a specific 
fiducial cosmic time (we can also identify Rq = H^^ when 
specifying the radial coordinate gauge through Ri). The 
initial conditions for this system are basically T-Lqi, 
and their radial gradients, which makes their specifica- 
tion very intuitive and practical, as these initial value 
functions could tend at a given asymptotic limit to the 
Hubble and Omega factors in a FLRW background at a 
fiducial cosmic time. The system (|117p can also be sup- 
plemented by the differential equations (|104p for L and 

r. 
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B. A dynamical system 



The system l|117p clearly suggests a dynamical systems 
approach, as discussed in 44] for FLRW and Bianchi 
models. This approach is based on the fluid flow equa- 
tions f32| - [33 | and an Omega parameter constructed as 
the ratio of density to the the squared expansion scalar 
li? = O'^/O. As shown in [i^l, if the evolution parameter 
is defined by d/dS^ — {l/'H)d/dT the Raychaudhuri equa- 
tion decouples from the remaining evolution equations, 
leading to a reduced system that can be analyzed qualita- 
tively: critical points, invariant subspaces, etc. For LTB 
models associated with evolution equations like (|117p . it 
is more convenient (see f29j) to define ^ in terms of H in 
(|118p by means of the coordinate transformation |29j 



T(e,r), 



r, 



^ = In H 



(119) 



so that for every scalar A(t, r) = A(t(^, r),r) — A(^,r) 
and derivatives associated with the 4- velocity flow d/dr 
(with r constant) become 



'dA 


OA 


-de 









H 



d_ 



)_d_ 
H97 



_ dA 
1 d 



Hq cdt 



(120) 



Hence, the three equations ()117b|) - (jll7dp become inde- 
pendent of H and effectively decouple from (|117ap . lead- 
ing to the reduced system 



dA 



(f2 - l) , (121a) 

(i + <5(«)) ,5(«) - i (^(«' + a) niuih) 

(121c) 



which is very similar to that analyzed in [29.] (where 
g^j^j dimensionless S were used instead of A and 
S^^''). Notice that this dynamical system is character- 
ized by a 3-dimensional phase space parametrized by 
{f2. A, (5('")}, which contains the FLRW dust case as an 
invariant set given by the line [Cl = n(^), A = 0, J^^-* — 
0]. Therefore, this dynamical systems approach also pro- 
vides a nice and appealing theoretical context to justify 
the role of Cl as an inhomogeneous generalization of the 
FLRW dust Omega parameter. 



XIV. REGULARITY OF SPECIAL LTB 
CONFIGURATIONS. 

A. Closed elliptic models and regular zeroes of 
R' = 0. 

In closed elliptic models the ^T[t] are homeomorphic to 
S'^ (see Appendix I A 3[) . There are two symmetry centers. 



at r = and r = Tc- Since i?(<, 0) = R{t,rc) = for 
all t, then i?'(rtv) — where < rtv < Vc must hold 
for all t. The regularity conditions (jA6P and (|A7|) imply 
that for all r the sign of R' is the same for all ■^T[t]. 
Hence, ([9]) requires the function Ri{r) to be selected so 
that Ri{0) = R,{rc) = and R'^in^) = 0, with 



i?- > and < J" < 1 for < r < rtv, 
i? ■ < and - 1 < J" < for rtv < r < 



(122) 



where J" = ±Vl + E, with J"(0) = 1, J"(rc) -1 and 
J-'(O) = P{rc) = -F(rtv) = 0. Since i?, R' and J" are 
bounded, if ^ holds then the maximal coordinate range 
is < r < Tc in all regular '^T[i]. 

Equations ()122|) provide the coordinate ranges where 
the "-h" or "— " signs hold in specifying the Hellaby-Lake 
conditions (jlip in the conventional variables. However, 
the ± signs are no longer needed with the initial value 
parametrization, as it is clear that fulfillment of (|72|) (the 
equivalent of (jlip ) implies by its construction that t'-^^^^ and 
t'^ji must have, respectively, the opposite and same sign 
(and a common zero) as R[/ Ri, and from (jA6p . this sign 
wiU be the same for R' / R at aU '^T[t]. 

If only the first sign equality in ((9]) holds, then there 
would exist a range of r for which 



sign(i?0 = sign(Af') ^ sign(J-), 



(123) 



so that M' and R[ have a common zero, but either F 
has a zero that is not common to that zero, or J- has 
no zeroes. In these situations, p in remains bounded, 
and so there is no curvature singularity, but we have a 
surface layer discontinuity at r = rtv [3, US E3] ■ This 
is the reason why when R' vanishes regularly at some 
r — rtv , this worldline must lie, either in an elliptic model 
or in an elliptic region. Besides the surface layer, (|123p 
also implies an ill-defined metric coefhcient in (P) and 
P9p (which leads in turn to an ill-defined proper radial 
length between different comoving dust layers becomes). 

Once these issues are taken under consideration, initial 
conditions complying with the Hellaby-Lake conditions 
([7^ can be provided for this class of models, whether for 
analytic/qualitive or numeric work. 



B. Simultaneous big— bang. 

In general, the initial curvature singularity given by 
(|58aP is not simultaneous, since it is marked by the curve 
ct = ctbb(r) in the (ct, r) coordinate plane. For parabolic 



(m) _ 



0, 



models, the condition ci^^ = in (|63ap implies 5, 
and thus m^j = 0, which corresponds to the FLRW limit. 
However, non-trivial hyperbolic and elliptic LTB models 
follow by settin g th e initial value functions so that ctbb is 
a constant (see [lO])- In order to examine these cases, we 
use ([5^ . ([57)1 and (|63bp to rewrite F in (pD|) in terms of 
ctbb and ct^i^ , which after setting t^^ = and tbb = ^1? = 



18 



constant, leads to 



r = 1 + 3((5,|"^ - jf')) - c{t - t\X') ( sy" - ^s. 



where Hq follows from psp and 



(0)^ I j^{,n) _ 3^(fc) 

(124) 



™^ (125) 



with Xt, yt, Zh and Ze given by (j49]), ^ and ([54l) . 
Since and yi depend only on ruqi and kqi, a simulta- 
neous big-bang t^^ = implies the following constraint 
on these initial value functions 



r.(i - / ^h{xt)/yt hyperbolic 

^ ' ' 1 Ze{xr)/y, elliptic 



(126) 



which, in turn, implies (by differentiation) the following 

constraint between the (5|™'' and 5,-'^'' (hence, between the 
gradients m^^ and /c^J 



1-^5^0(^,-4?) 



(127) 

As a consequence of p26p . a simultaneous big-bang al- 
lows us to prescribe only one of the two functions m,qi 
and kqi, the other one must be found by solving the al- 
gebraic constraint (|126p or solving (I127P as a differential 
equation for dkqi/dmqi (the function Ri{r) remains free 
and can be fixed as a radial coordinate gauge). Since 
we need both functions mqi and kqi to solve (jl02p or 
P03p . finding previously a numeric solution of the alge- 
braic constraint (|126p is not a problem (it is part of the 
work needed to provide appropriate initial conditions). 
However, assuming that iriqi is selected and bearing in 
mind that 



< Zh{xi) for Xi > 0, 

< ZeiXi) < TT for < Xi <2, 



(128a) 
(128b) 



it is evident that the hyperbolic branch of (|126l) does not 
appear to be restrictive at all to find \kqi\. For the elliptic 
branch, a necessary condition for a solution of (|126p that 
yields kqi is 



< C{U - il?)fc^r < ^TTT-qi^ 



(129) 



which, again, does not seem to be restrictive (see pC| for 
comparison). Still, since we cannot solve (|126l) analyti- 
cally, nor invert (|125p to find L, the case of a simultane- 
ous big-bang may be easier to handle analytically with 
the conventional variables M, E and R by means of ([^5]) 
or ([55]) . since we only need to fix t^,^ = t^^ and then pre- 
scribe the functions M and E taking care to fulfill the 
Hellaby-Lake conditions PT!]) and (fTT|) . Nevertheless, we 
can still use niqi and kqi to examine qualitatively these 
regularity conditions in these models through (|124l) . 



• Hyperbolic models 

We examine (|124p in the limits L w and L oo 
along constant but arbitrary r (so that the time 
dependence is concentrated on L). Consider the 
following series expansions of Hqcit — ti?): 



H5c(t-il?) 
for i < 1, 

HqC{t-t^^}) 

for L > 1. 



1/2 



2 X 
- + 

3 15 



L -1-0(^3/2) 



(130a) 



l + ^^i^ + 0(L-) 



XiL 



(130b) 



where we used (gH]), ^ and (fT25l) . Taking the 
leading terms in these expansions, F in (|124p be- 
comes 

F«l + (5|"^ (i«l), 

r«i + ^#-' (i»i), 

leading to the following conditions for the fulfill- 
ment of (ISTt 



which are the same as 
• Elliptic models 



,5("^ > -1. 



with t' = 0. 



(131) 



The behavior of F in around i « fori « t^°^ 
follows from an expansion similar to (ll30bP . but 
using Ze instead of Zh- This yields also Hqc{t — 

tl°}) w 2/3 -I- 0{L), and lead to F « 1 -I- S^'^l 
Considering now L sa but in the collapsing phase, 
so that t « i„„„, we have 'Hqc{t — t^?) —oo. 
Hence, F in this limit takes the form 



F«3|H,|c(i-ti?) (4"^-^<5f) 



(132) 



while F for t « takes the form . Therefore, 
the conditions for the fulfillment of ((6T|) are simply 



s-(m) 



> -1, 



(133) 



Notice, by comparing with ([75)) . that the second 
condition above is equivalent to ci'^j, > 0, hence 
(|133p are simply ([7^ with t^,,,^ = 0. It is important 
to notice that a simultaneous big-bang t^-^ = does 
not imply a simultaneous collapsing singularity or 
maximal expansion (t'^,j and t[_^^^^ are not zero). 



C. Simultaneous maximal expansion and 
simultaneous collapsing time. 

These elliptic FTB configurations follow by setting the 
maximal expansion and collapse times, i,„ax, ^con, in (|57p 
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to constants. They have been examined in [T3|, hence 
the reader can compare the treatment of these cases in 
this reference with our treatment in this subsection. 

The condition for t[^^^^ = in elhptic models follows 
directly from ([57)) as 

Ci„ax = Cti +70, TT - Ze{Xi) = JoUi, (134) 

where 70 is a positive constant and Xi, iji and are 
given by (f54| and (|56l) . As in the case t^^^ =0, we can only 
prescribe one of ruqi and kqi , with the other one obtained 
by solving numerically the constraint (jl34p . If niqi is 
prescribed, then a necessary condition for the existence of 
a solution for kqi follows directly from (|128b[) and (|134l) : 

< 7ofcg/^ < Trmgj, (135) 

which looks like ()129p and also does not appear to be 
restrictive at all. Still, it is necessary to work (|134|) nu- 
merically or to make further assumptions on rUqi or kqi to 
get more information. The locus of the big-bang and col- 
lapse singularities are now given by inserting (jl34p into 

m 

TT TT 

cibb = cti + jo , ct,„„ = cU + 70 H , (136) 

which indicates a time symmetric location of tbb and t^^^ 
with respect to t-^._^^. The gradients of tbb and ^^011 are 

_|=^.^(,«_|,<.). ,,37, 

where we used ([55]) and (PT|) . 

If i'^j, = 0, then ((57)) implies the following constraint 
similar to ()134p 

ct„n = cii + eo, 27r - Ze(a;i) = eoy^, (138) 

where eo is a positive constant. Again, we can only pre- 
scribe one of wLqi and kqi, and obtain the other by solving 
(|T38)) . From ()128bp we have tt < 27r-Ze(xi) < 27r, hence, 
if we prescribe niqi a necessary condition for finding kqi 
as a solution of (I138P is 

^3/2 

TT < £0^^ < 27r, (139) 

which, again, does not seem to be restrictive, though fur- 
ther information requires either numerical work on ()138p 
or making assumptions on rUqi or kqi . The big-bang and 
maximal expansion times and their gradients are 

27r 

cibb = cti + €0 , 

Vi 

TT 

ct,„^^ = cti + 60 , (140a) 

g^bb ^ ^ /^^A(fe) _ 



By comparing ([T37l) and (I140bl) with ([72]), then suffi- 
cient conditions to fulfill (|CT)) in either case ct^^^^ = and 
ct'caii — are simply ()133p . As with the case ct'^^^ = 0, 
it is not problematic to work out these cases in solving 
(|102l) or (|103l) numerically. As in the case ^^b = Oj these 
cases are easier to handle analytically with the conven- 
tional variables, since for a given choice of M and E, we 
have Hi = |i?|^/^/M, and thus cttb follows directly from 
(I136P or ()140ap . Of course, the three free parameters 
must comply with the Hellaby-Lake conditions (|TT)) . 



D. Mixed hyperbolic/elliptic configurations. 

LTB models admit more than one kinematic class in 
their full radial domain. This type of "mixed" configu- 
rations can be constructed either by smoothly matching 
LTB regions with one kinematic class to regions of an- 
other (as in [23]), or simply by choosing the initial value 
function kqi (which determines the kinematic class) so 
that it changes sign in its radial domain. 

A particularly interesting mixed configuration is an el- 
liptic region surrounded by a hyperbolic exterior given 
by the choice 

{> for < r < Tf,, elliptic region 
= for r = ri,, interface . (141) 

< for r > r^, hyperbolic region 

The fact that (5^'^'' — 00 if fcgi — (from its definition, 
see Appendix I A 4p . and thus (J''"'-' — )■ —00 for all t, signals 
a potential problem in using the quasi-local variables to 
study an elliptic/hyperbolic model. However, we notice 

that (JC") does not appear in (ITU^ nor in (ITU5)) . Also, (5,^^ 
only appears in the initial conditions (llOSp in the form 
kqi5f\ which docs not diverge as kqi — > 0. Hence, this 

behavior of 5^^'^ has no consequences for the numerical 
integration of these evolution equations for these config- 
urations. 

For analytic or qualitative work, the limit kqi — ;> must 
be handled carefully . While regularity conditions for 

these mixed models follow also from (|60)) . even if — > 
—00 as kqi — 0, we cannot simply set kqi = in ()60p . 
Instead, we examine F at constant values of r as r rj, 
(or kqi — > 0). The following series expansions around 
kqi — hold in this limit 

Hqc{t-ti) « ^(^l--L_^+0(kq,), (142a) 
1 - ^ « 1 - + 0(^9^)' (142b) 
so that F in (IMl) takes the form 

r^l+(^l-j^yt^+Oikq,), (143) 
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which (up to leading order) is the form of F in (l59l) . 
Hence, F passes smoothly from its elliptic to its hyper- 
bolic form, approaching the form for parabolic models 



k'{t, n) < and kq{t, n) = hold for all t: 



wherever k, 



for 



' n- 

Since the limit kqi — implies — > in either the 
elliptic (r < r^) or hyperbohc side (r < r^), we should 
also obtain for L the approximated parabolic form (|46)) 
for r ~ Tb from the analytic solutions ([50]) and ([55)) in 
the limit — )■ 0. Expanding Zh{u) and Ze{u) around 
M = yields 



Zh{u) 



,3/2 



^ ^ 3 6 



V2 



,5/2 



(144a) 
(144b) 



Hence, taking only the leading term in Zh and Ze and 
substituting into ([5D|) and the expanding phase of ([55]) 
yields for both elliptic and hyperbolic sides 

L"^'^ ~l + ^^^^^c{t-U) + 0{x^), (145) 



which up to the leading term coincides with (l46l) . We 
did not consider the collapsing phase in (|55|) because, 
from (l57l). we have cCx — ^ 



->■ oo as fcqi ~> 0, hence the 
worldline marking the interface r = r?, is contained in 
the expanding phase for all t. 

The conditions for the fulfillment of (|5T|) can be ob- 
tained jointly for the elliptic and hyperbolic regions, 
bearing in mind that for r « we have ^^''^ — >■ — oo but 
also F and L take the forms ^15)) and For both 

the elliptic and hyperbolic regions we have as r — >■ r;, 




(146) 



hence iJ-™^ < (with iJ,-™^ = only at r = 0) is a suf- 
ficient condition for ct'^^^ < in and (17^ to hold in 
the full radial range (note, from ([5T|) . that this condition 
is equivalent to m^j < 0). Since ct ^^uir) — >■ oo as r — > r?,, 
the collapsing singularity is contained entirely in the el- 
liptic region, hence the existence of the hyperbolic region 
for r > rb has no effect on the condition ci' > in (17^ . 

The only remaining condition is J^'^'' > —2/3 in the 
hyperbolic region. Evidently, this condition does not hold 



in the limit r — > rf, because S, 



(k) 



-oo, but in this limit 



F has the form (|143l) . and so (|61l) is not violated. As 
long as Jp-* remains negative we have /c^j > (from ([31]) 
with kqi < 0), and thus E' > holds, hence condition 
([T0| (equivalent to (|6T|) ) also holds. In the elliptic side, 
we have kqi — > with k'^^ < 0, which implies from (jB3p 
that i^' > and fc,; = ^72.^/6 < near rt,. In fact local 
spatial curvature k = is negative for all times near 

r = rfj, even if fc^ > all the way up to r = rb. This 
follows directly by applying the integral property (j27b[) 
with r = Tf, as integration limit and considering that 



k{t,rb) 



i?3(t,rb) 



k'{t,x)R^{t,x)dx <0. (147) 



This provides an example of how elliptic dynamics does 
not (necessarily) imply positive local curvature (see Ap- 
pendix B). 

Notice that if instead of a hyperbolic region in (|141l) . 
we have a parabolic region (fcg = 0) for r > rf,, then the 
results above still apply to the elliptic region in < r < 
r^. Hence, we would still have fc < near r = in this 
region, and thus fc < would also hold in the parabolic 
region. This is an example of how local spatial curvature 
can be negative in a parabolic region not containing a 
symmetry. 



XV. CONCLUSION AND FURTHER WORK. 

We have provided a comprehensive examination of 
LTB dust models in terms of quasi-local scalars and their 
fluctuations, which are covariant objects that can be re- 
lated to covariant scalars in the "fluid flow" or "1-1-3" 
formalism. The motivation, and contents of the article, 
together with a summary of previous work, concepts and 
ideas, have been given in detail in the introduction. As 
we have shown throughout the article, the initial value 
parametrization that emerges from these scalars is use- 
ful in analytic, qualitative and numerical studies of the 
models. The "fluid flow" evolution equations for these 
scalars are an appealing and practical alternative to the 
analytic solutions conventionally used in the literature, 
as they are fully general and at the same time techni- 
cally simple: they can be effectively handled as ordinary 
differential equations, and thus can be very handy for a 
numeric treatment of the models. These evolution equa- 
tions can also be understood in terms of a gauge invariant 
and covariant (GIC) perturbation formalism, consistent 
with the covariant fluid flow approach of Ellis et al (32| - 

tas well as the traditional gauge invariant approach 
. Under this formalism, the dynamics of LTB mod- 
els can be cast in terms of the dynamics of spherical, 
non-linear GIC perturbations on a FLRW background. 

Regarding analytic and qualitative work, the quasi- 
local scalar representation and its associated initial value 
description provide an appealing theoretical context to 
understand the conventional parameters of the models, 
leading to an understanding and appreciation of previ- 
ous work under a new perspective. This is specially evi- 
dent in the fact that all scalars can be expressed by sim- 
ple scaling laws that can be analyzed qualitatively, as 
well as in the formulation of regularity conditions (the 
Hellaby-Lake conditions 24:, 25, 39]) as restrictions on 
initial value functions. We have also looked at the prob- 
lem, proposed by Krasihski and Hellaby 0, |1| , of finding 
the existence of a unique LTB model that is consistent 
with the "mapping" of arbitrary radial profiles of den- 
sity or velocity at different cosmic times. Also, we have 
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provided a theoretical context for an Omega parameter 
for LTB models, which has been defined in the literature 
as an ansatz [l^, [l^ . Given the widespread use of LTB 
dust solutions for dealing with a wide variety of prob- 
lems and models in General Relativity and Cosmology, 
the new approach to these solutions presented here has 
a significant potential for applications. 

In this article we have mostly introduced, discussed 
and justified the use of a new scalar representation for 
studying LTB models. However, the formalism that we 
have presented is potentially useful in future work using 
these models, hence we are looking at its direct applica- 
tion in separate articles, on the make, dealing with new 
results on important theoretical and practical issues, such 
as: (i) the asymptotic properties and boundary condi- 
tions of scalars in the radial direction, (ii) the conditions 
for constructing LTB models with scalars having radial 
profiles of "clumps" or "voids", (iii) the conditions for 
the existence of a positive "back-reaction" in LTB mod- 
els, which would allow us to mimmic the effect of an 
accelerated cosmic ex pan sion in the context of Buchert's 
scalar averaging [21I, |22{ . Together with these articles 
under elaboration, we are also conducting a comprehen- 
sive numeric study of the models, based on the evolu- 
tion equations derived here. Finally, as shown by (45| . 
the extension of the formalism of quasi-local scalars can 
be suitably modified in order to apply it to the quasi- 
spherical Szekeres spacetimes. The main justification 
of the present article is then to motivate and induce 
researchers to consider the introduced variables in the 
study and application of LTB models, as well as to serve 
as the theoretical reference for the use of these variables. 



Appendix A: Regularity issues of the new variables. 

The regularity of LTB models in the conventional 
variables has been extensively discussed in the litera- 
ture 0, H, [H, m, [3^ . We examine this issue in terms of 
the quasi-local variables and their fiuctuations. 



Since £^(0) = M(0) = M'(0) = £^'(0) = 0, then mg,{0) = 
mi{0) and kqi{0) = ki{0) must hold, and thus 'Hqi{0) = 
HiiO). For any scalar A and its dual function Aq we have 

A{t,0) ^ Aq{t,0), A'it,0)^A'q{t,0)^0, (A2) 

The same conditions hold at a second symmetry cen- 
ter r — Tc in closed elliptic models (see Appendix I A 3p . 
Notice that a central singularity is now associated with 
L{ct,r) = 0. 



2. The Riemann tensor. 

Regardless of which parameters or variables we may 
use, the regularity at each spacetime point or surface 
in LTB models can be characterized by the continu- 
ity and finiteness of the Riemann tensor ^(o)(f)(g)(h) — 



e^j^-j , in an orthonornial tetrad basis e^^^ 



0, [131 • Considering the natural tetrad associated with 
the LTB metric in its form 



6(0) = " 



e(i) = 



-(2) 



1 

LRi sin 6* 



5° 

0^1 



(A3) 



the nonzero Riemann tensor basis components are given 
readily in terms of the new variables ruq and (J*-™' by 

'7^(o)(i)(o)(i) = 3m - 2mq = mq{l + 2,5^'^^), 
^(i)(2)(i)(2) = 3m - niq = mq{2 + 3(5(™)), 



7^ 



(0)(2)(0)(2) 



n 



(2)(3)(2)(3) 



2m„ 



(A4) 



These basis components also involve other scalars like 
n,Hq,k,kqimd 5^^\ which are related to m, m^ 
and 5^™) by the constraints ^ and the 

Hamiltonian constraint (|20p for the local covariant 
scalars 



1. Symmetry centers. 

We have only considered in this article LTB models 
having (at least) one symmetry center, which is a regular 
timelike worldline corresponding to a fixed point of the 
rotation group SO (3). This is a sufficient condition for 
integrals in (j24p to be finite in a domain 'd{r) defined by 
The symmetry center can be marked as r = (and 
r = Tc if there is a second one) . The following conditions 
hold: R{t, 0) = R{t, 0) = 0. Considering ^ and ^ 
and i?' — > 1 as r — > 0, we can mark the symmetry center 
by a zero of Ri (r) (distinct from the locus of a curvature 
singularity: L — Q) 25]. Hence 



2m - fc + = 2m - fc -I- {HqS'^'^^f 



(A5) 



L(rf,0)>0, r(ci,0) = l. 



(Al) 



Clearly, it is sufficient for the regularity of LTB models 
based on (|A4|1 that all the involved scalars (not only m 
and mq) are bounded and continuous, but it is not a nec- 
essary condition, since H and Hq could diverge because 
k and kq diverge, with m and m^ remaining bounded. 
Also, we note that the relative fluctuations S^^^^ might 
diverge in some cases when A and Aq are bounded. 



3. Topology of the space slices. 

Given the existence of (at least) one symmetry center, 
the admissible topologies (homeomorphic classes) of the 
^T[t] are 
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• "Open models" : ^T[t] homeomorphic to M.^. 
There is only one symmetry center, at r = 0. 
As a consequence of the regularity condition ([9]), 
R' > and E > —1 must hold for all r, and so this 
topology is compatible with regions or models of all 
kinematic classes (hyperbolic, parabolic or elliptic 
with —1<E<0). We discuss the asymptotic 
radial regime associated with this topology in a 
separate article. 

• "Closed models" : ■^T[t] homeomorphic to S'^. 
There are two symmetry centers, at r = and 
r = Tc. Since R{t,0) = R{t,rc) = for aU t, 
then there must exist a turning value of R so that 
i?'(rtv) = where < rtv < r^- 



Since ^ is valid in each '^T[t\ and F = ^/l + E does not 
depend on t, a zero of (characteristic of S"^ topology) 
or the fulfillment of > (characteristic of M'^ topology) 
will be common to all ^T{t). Hence, all ^T[t] belong to 
the same topological class and thus 



sign[i?'(t„r)] =sign[i?'(t,,r)]. 



(A6) 



must hold for every r and for any arbitrary pair of distinct 
hypersurfaces ^T[ti], ^T[tj]. As consequence, dH) can be 
given as an initial condition 



sign (i?-) — sign A/' — sign J^, 



(A7) 



specified on an arbitrary fiducial or "initial" hypersurface 



4. Possible blowing up of the S'^^K 

From their definition pip , it is evident that relative 
fiuctuations S'^^'^ will diverge under certain regular con- 
ditions {i.e continuous and finite A and Aq) if in the 
integration domain i?(r) of ([M]) Aq has a zero that is not 
a common same order zero of A. This situation occurs 
in various situations that do not violate regularity, such 
as the zero of Hq a,t t = t,„^^ as elliptic models or regions 
pass from expansion {Hq > 0) to collapse {Hq < 0), or in 
mixed elliptic/hyperbolic configurations where kq passes 
from positive to negative at some comoving radius r = rb 
(see section XIV D). However, the quantities J^'"^ J*^''-' 
and (5^^^ in all basis components in (IA4[) only appear in 
terms such as nin J^"), kqS'^''^ and nqS'-'^\ which do not 
diverge at a zero of rrig, kq or Hq. Hence, the blowing up 



of J*^^-* and S^'''> because of zeroes of kq or Hq does not 
affect the regularity of the solutions, as conveyed by the 
continuity and finiteness of (|A4|) . Nevertheless, in these 
cases it will be preferable to use AqS^^\ instead of S^"^\ 
for studying the time evolution of scalars (which justifies 
the use of (|103l) over (|102l) for numeric work involving 
expanding/collapsing regions or models). 



5. Restrictions of the radial range due to a 
curvature singularity. 

When we introduced the quasi-local scalars in (j24p , we 
assumed that the integration range 7?(r) defined by 
was fully regular. However, in general, the coordinate 
locus (|58ap of the central singularity (expanding or col- 
lapsing) is not simultaneous (i.e. not marked hy t — i^b = 
const., so that t[^^^ ^ 0). Consider the case of an expand- 
ing non-simultaneous singularity (big bang), marked by 
the curve [ti,-t,{r), r], with < 0, in the {ct, r) coordinate 
plane, so that L{ty,i,{r) , r) = (the collapsing singularity 
is analogous). The hypersurfaces '^T[t] for t < <bb are 
only regular for the semi open subset 



■d{r) = {a; I rbb < a; < r} C i?(r). 



(A8) 



where rbb is the intersection of tbb(^) and the constant t 
value associated with the '^T[t]. As a consequence, the 
integration range in the definition ([M)) for these '^T[i] 
must be t?(r), not i?(r). For a collapsing singularity in 
elliptic models we have exactly the same situation, but 
the involved hypersurfaces are those with t > icoii(^), 
with the lower radial bound given by r = rcoii marking 
the intersection oi t^^uir) and the ^T[i]. 

However, this range restriction has no consequence in 
the definition and usage of the quasi-local scalars Aq, 
as the involved integrals can be treated simply as stan- 
dard improper integrals. We define at each ^T[t] with 
t < t^bif) the incumbent integrals with their lower in- 
tegration limit as y = rbb + e, for an arbitrarily small 
e > 0, and then obtain the limit as e — > 0. Off course, 
since scalars like m, k, % diverge in this limit, m^, kq and 
T-Lq might diverge as well, but the functions are well de- 
fined and behaved in the range i){r). This restriction only 
prevents the Aq from taking values r < rbb for t < tbb('"), 
and so all results that involve these scalars can be trivially 
extended to include hypersurfaces '^T[t\ for t < t^b (r) and 
t > icon(r). 



Appendix B: Local vs. quasi— local spatial curvature. 

The kinematic class of LTB models is determined by 
the sign of the initial quasi-local spatial curvature kqi (or 
E = —kqiR^). From ([57)1 . the sign of kqi determines the 
sign of kq. It is evident that parabolic models (contain- 
ing a symmetry center) are spatially fiat, as kqi 



trivially implies fc„ 



k 



= h = o 

0. However, it is not obvi- 



ous if a given sign of kq in regular hyperbolic or elliptic 
LTB models also determines the sign of the local spatial 
curvature k = In order to examine this point, we 

rewrite (I41bp as 



kf, 



3r L2 



(Bl) 



We examine the relation between kq and k and their 
initial values for hyperbolic and elliptic models. 
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Negative spatial curvature: hyperbolic models 

In hyperbolic models or regions containing a center, 
r > and S^'^'^ > —2/3 hold everywhere (conditions ((6T|) 
and ([68])). therefore, as a consequence of (|37| and (jBl|) . 
we have: 



fc, < 



< 0, 



(B2) 



Hence, since all regular hyperbolic models or regions 
comply with fc^ < 0, then all these models or regions have 
also negative local spatial curvature. As explain further 
below, the converse is not true, as local spatial curvature 
can be negative in certain elliptic regions in which kg > 0. 

Elliptic models and positive spatial curvature. 

The relation between kq and k is more complicated 
in elliptic models or regions, since standard regularity 
(conditions (I5T|) and ([7^ ) do not place a lower bound 
on S'l'^K Hence, the possibility that S'l'^^ < —2/3 occurs 
cannot be ruled out, and so fc < can happen in regions 
where fc^ > 0. It is straightforward to show from ([5^ 

that ^l'^'' < —2/3 can only occur in a regular elliptic 
model or region (F > 0) if £" = for some < rtv < J' in 
a domain t?(r). As a consequence, we have (5f ^ > -2/3 
for all regular elliptic models in which E' < holds for all 
r (with E' = only at the symmetry center) , and so (|Bip 



implies in this case that kq > ^ k > everywhere. 

If there is a zero of E' at some r = r* , then 5^'''^ < —2/3 
will hold in some regions without violating regularity con- 
ditions. Since -1 < < and E{0) = 0, then for r « 
we must have E' < 0. Thus, the only possible configu- 
ration is: E' <0 and (jf^^ > -2/3 for < r < r*, with 

E' >0 and (jf ^ < -2/3 for r > r* . Further insight into 
this situation comes from rewriting (jl3p as 



E' = RR' {kq - 3fc) = -iRR' 



k'q/kq 



2R'/R 



(B3) 



If the regularity condition ([XT)) holds, then i?7Vl + E > 
0, and so for E'{r*) = to occur the local curvature 
k must decrease sufficiently to reach k{r*) = kq{r*)/3, 

hence E' > (or S'l'^^ < —2/3) leads to further decreas- 
ing of k, so that k{r*) < kq{r*)/3 holds for r > r* . A 
situation in which fc < occurs with kq > Q can easily be 
conceived: since k' and k'q are both monotonously nega- 
tive, we have a curvature clump and so kq > k, thus, if 
kq decays to zero sufficiently fast k might become neg- 
ative. This happens in the elliptic side of the mixed el- 
liptic/hyperbolic configuration examined in section XIV 
D. 
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We use the notation [t] to emphasize that t is a fixed 

arbitrary parameter. 

This equation corrects an error in f25|, where —2/3 < 
5'''' < was incorrectly assumed. 

The coordinate resemblance between FLRW dust cos- 
mologies and LTB models is further highlighted by using 
metric form (1391) for the latter. 



